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TECHNICAL  NOTE  3910 

THE  RESPONSE  OF  AN  AIRPLANE  TO  RANDOM 
ATMOSPHERIC  DISTURBANCES1 

o 

By  Franklin  W.  Diederich 


SUMMARY 


The  statistical  approach  to  the  gust-load  problem  which  consists  in 
considering  flight  through  turbulent  air  to  be  a  stationary  random  process 
is  extended  by  including  the  effect  of  lateral  variation  of  the  instan¬ 
taneous  gust  intensity  on  the  aerodynamic  forces.  The  forces  obtained  in 
this  manner  are  used  in  dynamic  analyses  of  rigid  and  flexible  airplanes 
free  to  move  vertically,  in  pitch,  and  in  roll.  The  effect  of  the  inter¬ 
action  of  longitudinal,  normal,  and  lateral  gusts  on  the  wing  stresses 
is  also  considered. 

The  method  of  analyzing  the  rigid-body  motions  is  similar  to  that 
used  for  analyses  of  the  dynamic  stability  of  airplanes,  in  that  the 
equations  of  motion  are  referred  to  stability  axes  and  expressed  in  terms 
of  conventional  stability  derivatives.  The  method  of  analyzing  the 
dynamic  effects  of  structural  flexibility  consists  in  an  extension  of  a 
numerical-integration  approach  to  the  static  aeroelastic  problem  and  is 
in  a  form  which  offers  the  possibility  of  calculating  divergence  and 
flutter  speed’s  with  relatively  little  additional  effort. 

The  mean-square  values,  correlation  functions,  and  power  spectra  of 
some  of  the  aerodynamic  forces  required  in  this  type  of  analysis  are  cal¬ 
culated  for  one  special  correlation  function  of  the  atmospheric  turbu¬ 
lence.  It  is  shown,  for  instance,  that  if  the  span  is  relatively  large 
• compared  with  the  integral  scale  of  turbulence,  the  mean-square  lift  and 
root  bending  moment  directly  due  to  the  gust  are  substantially  reduced 
when  the  differences  in  instantaneous  intensity  of  the  turbulence  along 
the  span  are  taken  into  account.  However,  if  the  motions  of  the  airplane 
are  taken  into  account  the  mean-square  root  bending  moment  may  be  increased 
as  a  result  of  these  differences.  Also,  the  mean-square  pitching  moment 
is  shown  to  be  substantially  increased  if  the  tail  length  is  relatively 

^This  report  represents,  except  for  some  minor  changes,  a  thesis 
submitted  in  May  195^  in  partial  fulfillment  of  the  requirements  for 
the  degree  of  Doctor  of  Philosophy  at  the  California  Institute  of 
Technology,  Pasadena,  Calif. 
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Now  at  the  Langley  Aeronautical  Laboratory,  Langley  Field,  Va. 
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large  compared  with  the  scale  of  turbulence.  Finally,  the  wing  stresses 
due  to  longitudinal,  normal,  and  lateral  gusts  are  shown  to  he  statis¬ 
tically  independent  under  certain  conditions. 


INTRODUCTION 


The  local  velocity  fluctuations  acting  on  an  airplane  flying  through 
atmospheric  turbulence  are  functions  of  time  defined  only  in  a  statistical 
sense  and,  hence,  constitute  a  stochastic  or  random  process.  Consequently, 
the  responses  of  the  airplane,  whether  they  are  motions  (linear  or  angular 
displacement,  velocities,  or  accelerations),  forces  (lift,  pitching  moment, 
bending  .moment,  and  so  on),  stresses,  or  any  other  phenomena  determined 
by  the  turbulence,  can  also  be  known  as  functions  of  time  in  only  a  sta¬ 
tistical  sense. 

This  paper  is  concerned  with  the  statistical  characteristics  of  those 
responses  which  have  a  bearing  on  the  loads  and  stresses  experienced  by 
the  airplane;  although  other  problems  such  as  those  relating  to  passenger 
comfort  or  to  the  stability  of  the  airplane  as  a  gun  platform  can  be 
treated  in  the  same  manner,  they  will  not  be  considered  here. 

The  first  approaches  to  the  gust-load  problem  which  use  the  statis¬ 
tical  techniques  developed  for  stationary  random  processes  appear  to  be 
those  of  references  1,  2,  and  3-  (An  earlier  investigation  concerned 
with  the  motions  of  an  airplane  in  turbulent  air  is  reported  in  ref.  4.) 

The  fundamentals  of  these  approaches  are  discussed  in  some  detail  in 
reference  1,  and  mention  is  made  therein  of  investigations  in  other  fields 
of  engineering  and  physics  that  have  dealt  with  the  problem  of  deducing 
the  statistical  characteristics  of  the  output  or  response  of  a  dynamic 
system  from  those  of  its  input.  The  mean-square  normal  acceleration  of 
a  rigid  airplane  free  to  move  in  one  degree  of  freedom,  namely,  vertical 
motion,  is  calculated  in  reference  2.  In  addition  to  being  rigid,  the 
airplane  is  implicitly  assumed  to  be  small  enough  for  all  its  components 
to  experience  the  same  gust  velocity  at  any  instant  of  time.  This  means 
that  the  span  of  the  airplane  must  be  small  compared  with  the  integral 
scale  of  atmospheric  turbulence,  which  on  the  basis  of  the  available 
knowledge  concerning  the  properties  of  the  atmosphere  (ref.  4,  for  instance) 
appears  to  be  in  the  order  of  several  hundred  to  1,000  or  2,000  feet; 
that  is,  the  span  of  the  airplane  must  be  less  than  about  100  feet. 

The  purpose  of  the  present  paper  is  to  extend  this  approach  to  large 
flexible  airplanes  free  to  move  in  all  directions.  As  used  herein,  the 
terms  "small"  and  "large"  airplane  refer  to  airplanes  which  are  very  small 
and  not  very  small,  respectively,  compared  with  the  integral  scale  of 
turbulence;  thus,  an  airplane  flying  in  a  wide  variety  of  atmospheric 
conditions  may  be  "small"  under  certain  conditions  and  "large"  under 
others.  Similarly,  the  terms  "rigid"  airplane  and  "flexible"  airplane 
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are  used  to  designate  airplanes  flying,  respectively,  at  speeds  far  below 
those  at  which  dynamic  and  aeroelastic  effects  become  important,  and  at 
speeds  at  which  these  effects  have  to  be  taken  into  account;  the  same 
airplane  can  thus  be  "rigid11  under  some  conditions  and  "flexible"  under 
others . 

Several  fundamental  assumptions  are  inherent  in  the  analysis  con¬ 
tained  in  this  paper.  In  the  first  place,  all  atmospheric  disturbances , 
motions,  and  structural  deformations  are  assumed  to  be  small  enough  to 
produce  forces  that  are  linear  and,  hence,  superposable.  Also,  the  turbu¬ 
lent  "input"  to  the  airplane  is  assumed  to  be  stationary  in  a  statistical 
sense;  that  is,  the  turbulence  in  the  plane  of  the  flight  path  is  homo¬ 
geneous.  For  the  large  airplane,  the  additional  assumption  is  made  that 
the  turbulence  is  axisymmetric  with  respect  to  vertical  axes,  a  condition 
less  severe  than  complete  isotropy.  The  statistical  characteristics  of 
the  turbulence  are  thus  assumed  to  be  invariant  under  a  translation  of 
the  space  origin  within  the  horizontal  plane  and  under  a  rotation  of  the 
coordinates  about  the  vertical  axis.  Finally,  Taylor's  hypothesis  to 
the  effect  that  time  displacements  are  equivalent  to  longitudinal  space 
displacements  is  assumed  to  be  valid. 

The  aerodynamic  forces  directly  due  to  atmospheric  turbulence,  which 
constitute  the  input  forces  for  the  dynamic  system  represented  by  the 
airplane,  are  calculated  in  the  first  part  of  this  paper  for  the  large 
airplane,  that  is,  for  the  case  where  the  spanwise  distribution  of  the 
intensity  of  turbulence  has  to  be  taken  into  account.  (The  effect  of 
spanwise  variation  of  gust  intensity  on  the  lift  has  been  treated  by  a 
slightly  different  method  in  ref.  5*) 

The  dynamics  of  the  rigid  airplane  are  considered  in  the  second  part. 
The  dynamic  system  is  now  represented  by  a  set  of  three  simultaneous 
ordinary  differential  equations,  rather  than  one  as  in  reference  2;  none¬ 
theless,  the  problem  of  calculating  the  required  transfer  functions  is 
still  one  of  simple  algebra. 

The  next  part  is  concerned  with  the  small  flexible  airplane  and  thus 
has  direct  application  to  fighter -type  airplanes  and  guided  missiles 
operating  at  relatively  high  speeds,  in  addition  to  serving  as  a  prelim¬ 
inary  to  the  last  part  of  the  paper.  The  dynamic  system  is  now  represented 
by  a  partial  differential  equation,  and  the  calculation  of  the  transfer 
functions  requires  the  solution  of  ordinary  differential  equations.  Once 
these  functions  are  calculated,  however,  the  statistical  techniques  are 
the  same  as  before,  as  a  result  of  the  fact  that  the  lateral  variation  in 
gust  intensity  is  ignored.  Either  modal  or  numerical-integration  approaches 
may  be  used  to  analyze  the  dynamics  of  a  swept-wing  airplane  with  arbitrary 
stiffness  and  mass  distribution.  Although  modal  approaches  have  usually 
been  preferred  in  the  past  for  similar  problems,  it  was  believed  that,  in 
view  of  the  highly  complex  nature  of  modern  aircraft  structures  and  the 
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advanced  type  of  computing  machinery  required  and  generally  available  for 
their  analysis,  the  numerical-integration  approach  would  be  preferable 
and  it  has,  therefore,  been  used. 

The  last  part  contains  the  analysis  of  the  large  flexible  airplane. 
The  statistical  problem  is  now  that  of  a  system  which  is  characterized 
by  a  partial  differential  equation  with  time  and  a  space  coordinate  as 
independent  variables  and  which  is  subjected  to  a  random  input  that 
varies  in  time  and  space,  so  that  more  is  required  than  the  transfer 
functions  from  the  gust  intensity  at  one  point  on  the  wing  to  the  stresses 
at  another.  The  particular  statistical  problem  presented  by  this  case 
is  considered  in  some  detail,  and  the  appropriate  transfer  fractions  are 
then  obtained  by  using  the  numerical-integration  approach  presented  in 
the  preceding  part  to  solve,  in  effect,  the  ordinary  differential  equa¬ 
tions  which  describe  the  wing  deformations  at  any  given  frequency. 


SYMBOIS 


A  aspect  ratio 

b  span 

C(k)  Theodorsen  function 

Cj^  lift  coefficient,  L/qS 

lift-curve  slope 

rolling -moment  coefficient,  L'/qSb 

C7  coefficient  of  damping  in  roll,  defined  as  positive  for 

^  positive  damping 

damping-in-pitch  derivative 

static  pitching  derivative 

c  chord,  parallel  to  plane  of  symmetry 

c  average  chord,  S/b 

c^  section  lift  coefficient  at  station  y,  Z/qc 


NACA  TN  3910 


5 


El  bending  stiffness 

e-|_  distance  from  section  aerodynamic  center  to  shear  center, 

fraction  of  chord 


distance  from  shear  center  to  section  center  of  gravity, 
fraction  of  chord 


% 

G(y,t)) 


distance  from  shear  center  to  the  midchord  point,  fraction 
of  chord 

distance  from  shear  center  to  the  3/4-chord  point,  fraction 
of  chord 

dimensionless  lift -inf luence  function  ( Green* s  function  for 
the  spanwise  lift  distribution) 


GJ 


torsional  stiffness 


g  acceleration  due  to  gravity 

H(o))  response  to  sinusoidal  oscillation,  Fourier  transform  of  h(t) 

h(t)  indicial-response  function 

Iv  mass  moment  of  inertia  about  X-axis 


Iy  mass  moment  of  inertia  about  Y-axis 

Jq,J^  Bessel  functions  of  the  first  kind,  order  0  and  1 

Kq,Kj_  modified  Bessel  functions  of  the  second  kind,  order  0  and  1 

Ki0  integral  of  Kq 

k  reduced  frequency,  <nc/2U 

k*  dimensionless  frequency,  aoL*/U 

L  lift 

L*  rolling  moment 

L*  integral  scale  of  turbulence 

l  distributed  lift  per  unit  distance  along  the  span 


pitching  moment 
bending  moment 

twisting  moment 

mass  (of  airplane,  unless  designated  otherwise  by  subscripts) 

distributed  mass  per  unit  distance  along  the  span 

distributed  twisting  moment  (about  axes  perpendicular  to  the 
plane  of  symmetry)  per  unit  distance  along  the  span 

dynamic  pressure 

radius  of  gyration  about  center  of  gravity;  longitudinal  dis¬ 
placement  corresponding  to  time  displacement  t 

wing  area 

T* 

scale  parameter, 
time 

mean  flying  speed 

longitudinal  component  of  gust  velocity 
horizontal  component  of  disturbed  motion 

lateral  component  of  gust  velocity 
weight  of  airplane 
vertical  component  of  gust  velocity 
vertical  component  of  disturbed  motion 

coordinate  along  mean  flight  path 

distance  from  intersection  of  elastic  axis  and  root  chord  to 
airplane  center  of  gravity 

tail  length,  distance  from  airplane  center  of  gravity  to 
aerodynamic  center  of  tail 
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modified  tail  length,  distance  from  intersection  of  elastic 
axis  and  root  chord  to  aerodynamic  center  of  tail 

coordinate  perpendicular  to  plane  of  symmetry 

coordinate  in  plane  of  symmetry  perpendicular  to  mean  flight 
path;  vertical  deflection 

inclination  of  chord  to  X-axis 

span  ratio,  b/L* 

autoconvolution  function  for  /(y) 

dimensionless  lift  distribution,  cc^jcC-^ 

dimensionless  lift  distribution  in  roll,  cc^/cC^ 

variable  of  integration  corresponding  to  y 
angle  of  pitch 
angle  of  sweepback 
mass  density  of  the  air 

time  displacement,  argument  of  time-correlation  function 
the  power  spectrum  $  in  the  case  of  axisymmetry 
two-dimensional  power  spectrum  (double  Fourier  transform  of  \jr) 
Sears  function  (unsteady- lift  function  for  gust  penetration) 
one -dimensional  or  point  power  spectrum 

two-dimensional  power  spectrum  (single  Fourier  transform  of  ijr) 
one -dimensional  or  point  correlation  function 
two-dimensional  correlation  function 
frequency  of  oscillation 
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Subscripts : 


e,e 

effective 

f 

fuselage 

r 

wing  root 

t 

tail 

u 

horizontal  component  of  turbulence 

w 

wing;  vertical  component  of  turbulence 

Matrix 

notation: 

t  ] 

square  or  rectangular  matrix 

n 

diagonal  matrix 

L  J 

row  matrix 

{} 

column  matrix 

Dots  over  symbols  indicate  derivatives  with  respect  to  time. 

AERODYMMIC  FORCES  RESULTING  DIRECTLY  FROM 
ATMOSPHERIC  TURBULENCE 


The  motions  of  a  rigid  airplane  depend  on  the  overall  forces  and 
moments,  whereas  the  stresses  of  a  rigid  airplane  and  the  motions  and 
stresses  of  a  flexible  airplane  depend  on  the  distribution  of  these 
forces,  as  well.  This  part  of  the  report  is  concerned  with  the  calcu¬ 
lation  of  the  integrated  and  distributed  forces  and  moments  directly 
due  to  atmospheric  turbulence  when  the  spanwise  variation  of  gust  inten¬ 
sity  has  to  be  taken  into  account.  (The  forces  and  moments  caused  by 
the  .motions  which  result  from  the  forces  treated  in  this  part  can  be 
calculated  by  conventional  methods  and  will  not  be  considered  here, 
although  the  combined  forces  will  be  considered  in  the  following  parts.) 
Thus,  this  part  serves  as  a  basis  for  all  the  material  presented  in  the 
later  parts  pertaining  to  the  large  airplane,  and,  hence,  the  fundamental 
notions  required  for  an  analysis  of  the  large  airplane  are  introduced 
here  and  discussed  in  some  detail. 

The  basic  approach  is  as  follows:  First,  the  instantaneous  value 
of  the  quantity  of  interest,  such  as  the  lift,  is  expressed  in  terms  of 
the  instantaneous  gust  intensity  at  a  point  and  a  suitable  influence 
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function.  For  the  lift  this  influence  function  can  and  -will  he  identified 
with  a  certain  lift  distribution  on  the  given  wing  in  reverse  flow;  the 
same  procedure  may  and,  in  the  case  of  the  rolling  moment,  will  be  fol¬ 
lowed  for  other  integrated  forces.  On  the  other  hand,  for  the  local  lift 
the  influence  function  is  the  Green's  function  for  the  three-dimensional 
unsteady-lift  problem  and  cannot  be  identified  with  an  easily  calculated 
lift  distribution  on  the  wing  in  reverse  flow.  Inasmuch  as  no  knowledge 
concerning  this  function  appears  to  be  available,  a  method  of  calculating 
an  approximate  Green's  function  for  this  problem  is  outlined  herein. 

The  required  influence  functions  for  integrated  effects  can  be  synthesized 
from  this  function,  and  if  the  associated  lift  distribution  in  reverse 
flow  cannot  be  calculated  conveniently  this  approach  may  be  preferable. 
This  technique  is  illustrated  here  by  means  of  the  bending  moment. 

The  next  step  consists  in  using  the  expression  for  the  instantaneous 
value  of  the  given  quantity  to  calculate  a  correlation  function  for  this 
quantity  in  terms  of  a  correlation  function  of  the  normal  component  of 
the  atmospheric  turbulence.  The  power  spectrum  for  the  given  quantity 
can  then  be  obtained  by  taking  the  Fourier  transform  of  its  correlation 
function.  This  power  spectrum  is  considered  herein  to  be  the  desired 
end  result,  because  the  mean-square  values  of  the  quantity  and  its  deriv¬ 
atives  can  be  obtained  from  it,  and  other  statistical  parameters  of 
interest  can  be  obtained  from  these  mean-square  values .  Several  alter¬ 
native  approaches  for  calculating  the  aforementioned  spectrum,  either 
from  the  correlation  function  or  directly  from  the  spectrum  of  atmospheric 
turbulence,  are  given  in  connection  with  the  lift  and  are  directly  appli¬ 
cable  to  other  quantities  as  well. 

In  this  part  of  the  report  the  assumption  is  made  that  the  influence 
functions  of  concern  can  be  written  as  products  of  a  function  of  time 
alone  and  a  function  of  distance  along  the  span  alone,  and  advantage  is 
taken  of  this  simplification  in  calculating  the  desired  spectra.  This 
restriction  is  removed  in  the  last  part  of  the  report,  where  the  influence 
functions  considered  cannot  be  separated  into  space -dependent  and  time- 
dependent  constituents,  and  so  the  approaches  outlined  there  are  general¬ 
izations  of  those  presented  in  this  part;  they  may  also  be  used  for  the 
calculation  of  the  spectra  of  the  quantities  considered  in  this  part  if 
the  assumption  concerning  the  influence  functions  is  not  valid. 


Definitions  of  Statistical  Parameters 

As  pointed  out  in  the  introduction,  the  intensity  of  the  vertical 
component  of  turbulence  w(t)  is  a  random  process,  so  that  the  resulting 
forces  can  also  be  known  only  in  a  statistical  sense.  The  purpose  of 
this  part  is  to  calculate  certain  statistical  properties  of  these  forces 
namely,  their  mean-square  values,  their  correlation  functions,  and  thei? 
power  spectra.  The  fundamental  principles  involved  in  statistical 
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analyses  of  the  type  considered  herein  are  expounded  in  some  detail,  and 
citations  of  the  literature  on  the  subject  are  given  in  references  1 
and  2.  These  fundamentals  "will  therefore  not  be  repeated  here.  However, 
both  for  the  sake  of  ready  reference  and  inasmuch  as  the  statistical 
terms  are  not  always  defined  in  the  same  manner,  the  forms  that  are  used 
herein  are  indicated  in  the  succeeding  paragraphs. 

The  time  average  of  a  time -dependent  quantity  is  designated  by  a  bar 
placed  over  the  symbol,1  and  is  defined  as  follows: 

PT 

£  =  lim  -km  /  !(t)  at 

T— >oo  J  -T 

The  assumption  will  always  be  made  that  this  limit  exists  and  is  invarianl 
under  a  translation  of  the  origin  of  time.  This  assumption  implies  that 
the  processes  considered  here  are  stationary  in  a  statistical  sense. 

The  mean  of  a  random  process  f(t)  is  defined  as  its  time  average, 
and  is  always  assumed  to  be  zero.  In  dealing  with  processes  with  nonzero 
mean  this  analysis  is  thus  pertinent  only  to  the  process  which  consists 
of  the  difference  between  the  original  process  and  its  mean  value.  Simi¬ 
larly,  the  mean-square  value  of  a  random  process  f(t)  is  defined  as 
the  time  average  of  the  square  of  the  process,  so  that 

—  PT 

f2  h  lim  1  f2(t)  dt 

T— >00  2T  J  -  T 

The  time-correlation  function  of  f(t)  is  defined  as 


+f(r)  =  f(t)  f (t+r ) 


so  that 


=  tf(o) 


^When  no  possibility  of  confusion  exists,  a  bar  is  also  used  to 
designate  a  space  average,  as  in  the  case  of  c  and  y.  Also,  for  the 
components  of  turbulence  the  correlation  functions  depend  primarily  on 
space  displacements  and  can  be  defined  by  space  averages.  For  the  sake 
of  consistency,  however,  they  are  considered  to  be  defined  by  time 
averages . 
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and  the  power  spectrum  of  f(t)  is  defined  as  the  Fourier  transform  of 
the  correlation  function: 

00 

9f(a>)  =  i  j  e~X(m  \|rf(T)  dT 

d  —00 


dr 


The  second  form  is  the  more  convenient  one  when  \|r^  represents  more 

nearly  a  space  correlation  than  a  time  correlation,  so  that  it  depends 
directly  on  a  space  displacement  r  =  Ut  rather  than  on  a  time  displace¬ 
ment  t.  For  the  purposes  of  the  analysis  presented  herein,  the  gust 
or  input  correlation  functions  have  this  property,  so  that,  for  instance, 
may  he  defined  as 


fw(r)  =  w(t)  w(t +£) 

By  virtue  of  the  reciprocal  properties  of  Fourier  transforms  and 
the  symmetry  of  ^(t),  which  is  a  direct  consequence  of  the  assumed 

stationarity  of  f(t),  f^  may  he  expressed  in  terms  of  the  spectrum  of 
f  as 


f2=J  ^ 


(1) 


If  f(t)  is  the  input  of  a  linear  system,  the  power  spectrum  of 
the  response  x(t)  of  the  system  is  related  to  the  power  spectrum  of 
f(t)  hy 


(2) 


where  H(oo)  is  the  transfer  function  of  the  system,  that  is,  the  complex 
amplitude  of  the  response  of  the  system  to  unit  sinusoidal  input.  Hence, 
H(cu)  is  also  the  Fourier  transform  of  the  indicial  response  h(t)  of 
the  system,  which  is  defined  herein  as  the  response  of  the  system  to  a 
unit  impulsive  input: 


00 


e-icut 


H(o))  = 


—00 


h(t)  dt 


(3) 
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where  the  lower  limit  could  he  taken  as  zero,  since  h(t)  is  zero  for 
t  <  0.  Conversely,  h(t)  can  he  obtained  from  H(ao)  by  means  of  the 
inverse  of  equation  (3): 


h(t)  =  —  e***^  H(cd)  da) 

J  _00 


(3a) 


The  .mean-square  value  of  the  response  can  then  he  obtained  by  inte¬ 
grating  its  spectrum.  Similarly,  the  mean-square  values  of  the  nth  time 
derivative  of  the  response  can  be  obtained  from  the  (2n)th  moment  of  the 
output  spectrum.  For  instance, 


(S)  ■  f0 

From  the  mean-square  values  of  these  derivatives  other  statistical  quan¬ 
tities  of  interest,  such  as  the  expected  number  of  peaks  of  the  response 
per  unit  time,  can  then  be  calculated. 


Lift-Influence  Functions  in  Unsteady  Flow 

At  any  time  the  lift  on  a  wing  which  results  directly  from  atmos¬ 
pheric  disturbances  can  be  expressed  for  an  unswept  wing  as 


where  h(t,y)  dy  is  a  lift-influence  function  which  represents  the  lift 
caused  by  an  impulsive  vertical  gust  of  width  dy  which  at  time  t  =  0 
impinges  on  the  wing  at  station  y. 

The  influence  functions  required  in  equation  (k)  are  difficult  to 
calculate  directly;  methods  for  obtaining  lift  distributions  on  wings  of 
finite  span  in  unsteady  flow  usually  require  numerical  solutions  which 
do  not  lend  themselves  readily  to  the  analysis  of  angle-of -attack  distri¬ 
butions  represented  by  delta  functions.  However,  by  virtue  of  the  reci¬ 
procity  theorems  of  linearized  lifting -surf ace  theory  (ref.  6,  for 
instance)  the  lift  influence  function  for  a  twisted  wing  in  indicial 
motion  is  equal  to  the  lift  distribution  on  that  wing  during  indicial 
motion  in  the  reverse  direction  with  uniform  unit  angle  of  attack.  The 
lift  distribution  in  indicial  motion  with  uniform  angle  of  attack  can 
be  calculated  relatively  easily. 


NACA  TN  3910 


13 


For  the  few  cases  for  which  calculations  have  heen  made  (namely, 
some  unswept  wings ),  this  lift  distribution  tends  to  be  substantially 
invariant  in  time,  except  for  overall  magnitude.  For  instance,  the 
calculations  of  reference  7  indicate  that  the  lift  distribution  of  an 
oscillating  rectangular  or  elliptic  wing  in  incompressible  flow  is  sub¬ 
stantially  independent  of  frequency,  so  that  in  indicial  motion  it  is 
substantially  independent  of  time.  This  simplification  may  not  be  valid 
for  swept  wings. 

The  lift  influence  function  can  then  be  written  as 

h(t,y)  =  i  hL(t)  y(y)  (5) 

where  y(y)  defines  the  steady-state  lift  distribution  for  uniform  unit 
angle  of  attack: 


7(y)  = 


and  where  h^(t)  describes  the  variation  of  the  overall  magnitude  of 
the  lift  as  a  function  of  time  after  entry  into  a  sharp-edge  gust  and 
may  be  written  as 


hL(t) 


CIuqS  d  k  Alt  \ 

~a~  «  2w^J 


In  turn,  kg  is  the  lift  response  to  a  unit  sharp-edge  gust  normalized 

to  a  steady-state  value  of  unity  and  is,  as  in  the  preceding  equation, 

usually  expressed  in  terms  of  semichord  lengths  traveled  Hi-.  The  time 

c/2 

derivative  is  taken  here  hecause  the  response  wanted  is  the  one  to  a 
unit  impulsive  gust  rather  than  a  sharp-edge  gust. 

The  Fourier  transform  HL(o>)  of  this  function  hL(t)  is  proportional 

to  a  function  $(k),  which  may  be  termed  the  generalized  Sears  function 
because  for  two-dimensional  incompressible  flow  it  is  the  Sears  function: 

CT  qS 

HLU)  =  . .  0(k)  (6) 


where 


J 


k  s 


(DC 

2U 


Ik 
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The  function  Hj/oi)  represents  the  complex  amplitude  of  the  lift  due 

to  sinusoidal  gusts  of  unit  amplitude  and  is  thus  the  transfer  function 
from  the  gust  to  the  lift. 

Actually,  only  the  absolute  square  of  $(k)  will  be  required.  For 
two-dimensional  incompressible  flow  the  following  approximation  is  given 
in  reference  1: 


1 

1  +  2itk 


(7) 


This  expression  has  the  advantage  of  simplicity,  although  it  is  somewhat 
in  error  compared  with  the  abolute  square  of  the  Sears  function  at  very 
low  frequencies,  a  fact  which  could  be  remedied  by  using  the  approximation 


1  +  ak _ 

1  +  ak(l  +  2itk) 


where  a  is  about  15  for  a  good  overall  fit  to  the  exact  expression. 
However,  the  behavior  of  the  Sears  function  itself  at  very  high  frequencies 
is  unrealistic,  because  its  absolute  square  goes  to  zero  as  l/k,  whereas 
for  any  nonzero  Mach  number  and  any  finite  span  the  absolute  square  of 
the  generalized  Sears  function  can  be  shown  to  tend  to  zero  at  least  as 

rapidly  as  l/k2.  As  a  result  of  these  discrepancies,  the  approximation 
given  by  equation  (7)  and,  for  the  same  reason,  the  absolute  square  of 
the  exact  Sears  function,  cannot  be  used  to  obtain  moments  of  the  lift 
spectrum,  that  is,  values  of  the  mean-square  derivatives  of  the  lift, 
although  they  may  be  adequate  for  calculating  the  mean-square  value  of 
the  lift  itself  in  many  cases,  particularly  when  absolute  accuracy  is 

not  required.  Whenever  possible  the  values  of  |0(k)|2  given  in  refer¬ 
ence  8  for  the  plan  form  and  Mach  number  of  interest  should  be  used. 


Mean-Square  Lift  and  Its  Spectral  Resolution  for 
the  Unsvept  Wing 

Basic  equations.-  The  correlation  function  of  the  lift  can  be 
expressed,  by  virtue  of  equation  (4-),  as 


A00  A00  nb/  2  pb/2  - 

=  J  co  J  »  J  b/2J  b/2  h(tl,yl)  h(t2^y2)  w(U(t-tl)’yl)  w(U(t+T-t2)>y2)  <^1  ^2  dtl  dt2  (8) 


3N 
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where  the  averaged  product  on  the  right  side  represents  a  velocity  corre¬ 
lation  function.  This  function  depends  in  general  on  both  space  and  time 
displacements.  However,  if  Taylor's  hypothesis  is  made,  the  time  dis¬ 
placements  are  equivalent  to  longitudinal  space  displacements.  The  veloc¬ 
ity  correlation  functions  are  then  functions  only  of  longitudinal  and 
lateral  space  displacements.  Thus,  for  homogeneous  turbulence, 

w(x+Ut,y)  w(x+|+U(t+r),y+T])  =  tw(|+Ur,rj)  (9) 

In  addition  to  Taylor's  hypothesis  and  the  assumption  of  homogeneity, 
the  turbulence  is  assumed  to  be  axisymmetric  with  respect  to  vertical 

axes,  so  that  is  a  function  only  of  +  t]2.  This  function 

is  the  ordinary  space-correlation  function  ilr  ,  so  that 

W 

tw(UT+|,T))  =  lrw^(UT+i  )2+712) 

and,,  hence , 


If  the  assumption  implicit  in  equation  (5)  is  now  made,  the  preceding 
equation  can  he  written  as 


where,  in  turn,  r(rj)  is  an  autoconvolution  of  y(y)  defined  by 


(12) 


T(r|)  =  - 


D  ~ 

2"^ 


y(y)  /(y+Ti)  ay 


v 
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The  validity  of  equation  (ll)  can  be  demonstrated  by  performing  the 
integration  in  the  y-^y^  plane  as  indicated  in  the  following  sketch: 


That  is,  integration  is  performed  first  over  y-j_,  with  a  variable 
rj  =  y^  -  y-j_  held  constant.  For  this  integration  is  constant,  so 

that  only  and  7(yp+ Tl)  are  involved,  and  the  result  is  a 

function  of  T]  which  is  one-half  of  the  function  r( rj )  defined  by 
equation  (12).  The  second  integration  is  then  performed  over  r\, 
yielding  equation  (ll)  except  for  a  factor  of  1/2.  In  this  process, 
only  the  part  of  the  square  above  the  line  y ^  -  y^_  is  covered.  How¬ 
ever,  by  a  similar  process,  the  part  of  the  integral  corresponding  to 
the  part  of  the  plane  below  this  line  can  be  evaluated  and  shown  to  be 
equal  to  the  first  part,  so  that  by  defining  T  as  in  equation  (l2), 
both  parts  are  taken  into  account  simultaneously  in  equation  (ll). 

The  quantity 

we2  = 

=  bj0  dT1 

may  be  considered  an  averaged  mean-square  vertical  component  of  turbu¬ 
lence;  \|rw  (Ut)  is  then  the  corresponding  correlation  function,  and  the 

Fourier  transform  of  the  latter, 
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V 


d(UT) 


(13) 


is  the  corresponding  power  spectrum. 

Once  cpWe (oo)  has  been  obtained,  the  power  spectrum  of  the  lift 

can  be  obtained  by  taking  the  Fourier  transforms  of  both  sides  of  equa¬ 
tion  (lO).  The  result  is 

cpL(o>)  =  |hl(oo)  I  cpWe(a>)  (110 


where  is  the  transfer  function  defined  in  the  preceding  section. 

The  mean  square  of  the  lift  can  then  be  obtained  by  integrating  its 
spectrum,  as  indicated  in  equation  (l). 

Equation  ( lb )  has  the  same  form  as  the  corresponding  equation  for 
the  case  where  spanwise  averaging  of  the  effects  of  turbulence  is  not 
taken  into  account  (see  eq.  (2)  and  ref.  l),  except  that  cpv(a>)  is  now 

replaced  by  v  (u>) .  Thus,  the  spectrum  of  the  averaged  turbulence  must 

approach  that  of  the  unaveraged  turbulence  when  the  span  approaches  zero, 
as  may  be  seen  to  be  the  case  from  equation  (ll)  and  the  definitions  of 
7(y)  and  I^t]). 

Two  alternative  approaches.-  The  defining  relations  for  cpWe(co), 

equations  (ll)  and  (13),  do  not  necessarily  represent  the  best  method 
of  calculating  it  in  any  given  case.  A  slightly  different  expression 
appears  to  be  more  convenient  in  general.  It  consists  in  substituting 
equation  (ll)  into  equation  (l3)  and  inverting  the  order  of  integration 
to  yield 


r(q)  ?w(cu,T))  dT] 


(15) 


where 


9w(o),ii) 


(16) 


9w(ao,0)  =  cpw(o>) 


so  that 
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A  third  approach,  which  has  certain  advantages  over  the  others,  is 
similar  to  the  one  that  has  been  used  in  reference  5*  In  this  approach 
the  assumption  of  axisymmetry  is  not  made  initially,  and  use  is  made  of 
the  spectrum 


3>w(A1'A2)  "  *2 


1  /■"  r 

OO  vj  —00 


d|  dn 


(17a) 


The  correlation  function  (Ut)  can  then  be  written  as 


f'w2  w  w  ^  iLL  «1[Xi0,W2(!r^  mxi’i2)  “1 


Substitution  of  this  expression  into  equation  (13)  and  then  interchanging 
the  order  of  integration  yields: 


pOO 

3S,(g<^)|?(^)|2  d*2 


where 


rb/2 

r(A)  =  i  /  7(7)  e_iAy  dy 

b  J- b/2 

and  is,  as  a  result  of  the  symmetry  of  y(y),  real  and  symmetric  in  A. 


Now,  if  the  turbulence  is  axisymmetric,  ®v  depends  only  on 
A  s  |a-|_2  +  A22,  that  is, 


(17b) 


so  that  the  expression  for  cpw  (cu)  becomes 


<P, 


we(a))  =ul  Mtu)2^2 rM 


dA 


(18) 
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where  $V(A)  is  related  to  ^(r)  by 

poo 

=  |  JQ  r  J0(Ar)  1V(r)  dr 

Although  the  spectra  used  in  the  last  two  approaches  have  been 
defined  by  expressions  involving  the  point  correlation  function  \|r  (r), 

they  can  be  expressed  equally  well  in  terms  of  the  point  spectrum  cp  (od), 

w 

so  that  if,  say,  an  experimentally  obtained  point  spectrum  is  to  be  used, 
it  need  not  be  transformed  into  a  correlation  function  before  it  can  be 
used  in  these  calculations.  The  required  relations  are 


and 


-  - 


CO 


a2u2 


,3/2 


dCD 


(20) 


where  the  notation  I  is  used  to  specify  that  the  finite  part  of 

the  integral  is  to  be  taken,  an  operation  which  .may  be  performed  by 
integrating  by  parts  and  ignoring  the  infinite  part,  so  that,  in  terms 
of  a  proper  integral. 


*  (-K)  -  _  suf  r  dCP^(cu)  da) 

V  *  4l  *■> 

Also,  the  function  can  be  obtained  from  cpw(o>,T))  by  means  of 

either  of  the  relations 

p,  CO 

=  rfjj  e~±M  dT) 

P  00 

=  Jr  /  dTi 


0 
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Results  of  calculations.-  In  order  to  illustrate  the  magnitude  of 
the  effects  under  consideration,  calculations  have  been  madfe  for  a  uni¬ 
form.  loading  /(y)  =  1  and  a  point  correlation  function  which  has  been 
used  in  references  1  and  2  and  appears  to  fit  experimental  data  (ref.  4) 
fairly  well  over  a  large  portion  of  the  significant  frequency  range, 
namely, 

Vr)  =  w2(l  “  ^|A)  e",r,/L* 


so  that 


%>)  = 


w^L* 

itU 


1  +  3k'2 

(l  +  k'2)2 


where  L* 


is  the  integral  scale  of  turbulence,  which  is  here  defined  as 


if  (r)  dr 


and  where 


This  correlation  function  has  the  drawback  that  the  .moments  of  the 
spectrum  associated  with  it  are  infinite,  so  that  it  implies  a  process 
with  infinite  mean-square  derivatives,  but  it  is  quite  useful  if  only 
the  .mean-square  value  of  the  process  itself  is  of  interest. 


For  uniform  loading , 
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and  for  the  given  correlation  function, 

'%;K  if 

lL 


?w(a),y)  = 


2t* 
w  L 


*U 


£  7  >  KJ^tk 

(l  +  k'  ) 


y 


-  ($’  777 


vhere  Kq  and  are  .modified  Bessel  functions  of  the  second  kind,  and 


_  jwV2 


,2 


Jtu2  /  t2\5/2 

(1  +  k* d) 


The  mean-square  average  gust  intensity  for  this  case  is  given  by 

"2  ~2  1  -  e~P 


we  =  w 


P 


where  p  =  b/L* ,  and  is  shown  in  figure  l(a).  The  correlation  function 
tw  is  given  by 


-  f 


a  sinli-1  Kq^uj  sinh-i  +  — ^e-cr  - 


where  ff  =  r/L*,  and  where  Kq(c>0)  and  K^(crj9)  are  incomplete  modi¬ 
fied  Bessel  functions  of  the  second  kind  defined  by 

Kn(ff>e)  H  r  cosh  n3  e“a  cosh  «  d3 
JG 

The  spectrum  <Pw-e(m)  is  given  by 

%(“)  =  ^  §f  - - - 3  <  3k'%]jlTk^  [^(p^k'2")  - 

P  (l  +  k*2)  1 


pyi  +  k’2  ^(pjl+k'2^  +  (l  -  3k’2) 
p2(l  +  k'2)  K0(pjl+k'2) 


2  -  2p||l  +  k'2  K^pjl+k'2) 
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where  KlQ(x)  is  the  integral  of  Kq(x): 

Ki°(x)  =Io  Mxi)  dXl 

The  functions  and  epw  (normalized  with  the  averaged  mean- 

e  e 

square  turbulent  velocity)  are  shown  in  figures  2  and  3,  respectively. 

The  effect  of  the  span  ratio  on  the  normalized  correlation  function  of 
the  averaged  turbulence  may  be  seen  to  be  relatively  small,  so  that  the 
effect  on  the  unnormalized  correlation  function  is  primarily  the  decrease 

in  overall  level  given  by  the  ratio  of  we2  to  wl  A  similar  statement 

may  be  made  for  the  power  spectrum.  If  the  power  spectrum  were  not 
normalized  the  averaging  effect  of  the  span  would  tend  to  reduce  the 
intensity  of  the  spectrum  at  all  frequencies,  but  the  high  frequencies 
would  be  attenuated  much  more  than  the  low  ones,  as  might  be  expected. 

In  fact,  although  the  unaveraged  spectrum  decreases  as  oT2  at  high 
frequencies,  the  averaged  spectrum  decreases  as  of"* . 

The  asymptotic  values  for  b/L*— ><»  are  shown  in  figures  2  and  3 
in  order  to  indicate  the  nature  of  the  functions  considered  here  when 
the  scale  of  turbulence  is  small  compared  with  the  span,  as  may  be  the 
case  for  a  wind-tunnel  model  responding  to  natural  or  artificial  tunnel 
turbulence,  or  for  a  buffeting  wing  or  tail  surface,  although  this  con¬ 
dition  is  not  of  practical  concern  for  the  gust-load  problem. 

The  power  spectrum  of  the  lift  is  equal  to  the  product  of  the  power 
spectrum  *we(  cjq)  and  the  absolute  square  of  the  transfer  function  H^(m), 

as  indicated  in  equation  (l4) .  Inasmuch  as  this  lift  is  not  an  end  in 
itself  but  only  one  of  the  parameters  that  enter  into  the  calculations 
of  the  motion  of  the  airplane,  its  mean-square  intensity  is  of  little 
practical  significance;  its  spectrum  is  the  quantity  needed  in  further 
calculations.  However,  if  the  mean-square  intensity  is  wanted  for  any 
reason  it  can  be  obtained  by  integrating  the  spectrum.  Thus,  for  instance, 

the  approximate  expression  for  |$(k)|2  given  by  equation  (7)  and  the 
spectrum  9W(<^)  used  for  the  preceding  calculation  yields  the  mean-square 

lift: 
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where 


s 


and,  in  view  of  the  observation  that  much  of  the  turbulent  energy  is 
contained  in  a  region  for  which  the  span  has  a  very  small  effect  on  the 
(normalized)  spectrum,  this  equation  should  serve  to  furnish  an  approxi¬ 
mation  to  the  mean-square  lift  for  nonvanishing  span,  provided  we2  is 

used  instead  of  w^.  However,  the  mean-square  values  of  the  derivatives 
of  the  lift  cannot  be  calculated  in  this  simple  manner,  because  the 
deviation  of  the  normalized  averaged  spectrum  from  the  point  spectrum  at 
high  frequencies  cannot  be  ignored  in  calculating  the  moments  of  the 
lift  spectrum. 


Mean-Square  Lift  of  the  Swept  Wing 

For  the  yawed  or  sideslipping  unswept  wing,  equations  (10)  and  (l4) 
for  the  lift-correlation  function  and  spectrum  are  still  valid  if  an 
appropriate  lift-influence  function  is  used,  and  if  the  correlation 
function  for  the  averaged  turbulence  is  defined  by 


V(ut)  =  r — ~ — r 

e  b  cos  A 


■1)  cos  A 


r(rj)  tw(j(UT+r,  tan  A)2+t)2^  dT] 


where  r( rj )  now  pertains  to  a  lift-distribution  function  ?(y)  which  is 

appropriate  for  yawed  motion  and  is  defined  for  -  ^  cos  A  =  y  =  ~  cos  A. 

The  mean-square  averaged  intensity  of  the  vertical  component  of  turbulence 
is  then 


1)  cos  A 


b  cos  A 


rU)  * J— 2— ) 

\cos  A/ 


dr] 


1 

b 


cos  A)  \|/w(t)'  )  dll 


T 


where  V  =  — 1— .  Thus,  this  mean-square  intensity  is  unaffected  by  the 
cos  A 

yawing  process,  except  for  the  slight  change  which  results  from  the  change 
in  r( rj ) ,  although  the  spectral  resolution  of  the  averaged  intensity 
changes  in  the  process. 
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For  the  swept  wing  both  y2  -  y1  and  |  y2  |  -  | y1 1  occur  in  the 
integral,  so  that  the  reduction  of  the  double  integral  for  \|r  (Ut)  to 

a  single  integral  (see  eq.  (ll))  cannot  be  affected  so  simply.  The 
double  integral  for  the  swept  wing  is 

+We(UT)  =  ^  J_t/2  7(yL)  7(y2)  +vypT+(|y2|-|y1|  )t«m  A]2+(y2-yi)2  j  dyi  dy2  (21) 


From  this  integral,  by  using  rectangular  lift  distributions  and  the 
aforementioned  point  correlation  function,  the  mean  square  of  the  averaged 
turbulence  as  well  as  the  corresponding  correlation  function  and  power 
spectrum  have  been  calculated  by  numerical  integration  for  various  sweep 

angles  A,  the  ratio  - being  maintained  at  0.5*  (The  decision  to 

COS  ' 


rather  than  p  constant  was  reached  as  a  result  of  the 


hold  — — 
cos  A 

foregoing  analysis  of  the  yawed  unswept  wing,  which  indicated  that  the 
effects  of  sweep  should  be  minimized  in  this  manner.)  The  results  for 


we2  are  shown  in  figure  1,  and  the  effect  of  sweep  on  we2  is  seen  to 
be  small  for  this  comparison.  The  calculated  correlation  functions  and 

spectra  (normalized  with  respect  to  ve2)  are  not  shown  because  they 

agreed  with  those  for  A  =  0  within  less  than  1  percent  for  most  values 
of  Ut  and  k!,  respectively. 


Mean-Square  Rolling  Moment 

In  the  preceding  sections  the  averaging  effect  of  the  span  has  been 
shown  to  consist,  essentially,  in  reducing  the  effective  intensity  of  the 
turbulence  sensed  by  the  wing;  thus,  it  only  modifies  the  forces  present 
on  a  wing  of  small  span.  If  the  analysis  is  extended  to  the  rolling 
moment,  however,  a  new  phenomenon  appears.  When  a  wing  is  so  small  relative 
to  the  scale  of  turbulence  that  at  any  instant  all  of  its  points  experience 
the  same  turbulent  velocity,  the  wing  experiences  no  rolling  moment  as  the 
result  of  the  direct  action  of  turbulence  (although  it  may  experience  a 
rolling  moment  indirectly  as  a  result  of  the  rolling  and  yawing  motion 
caused  by  the  lateral  component  of  the  turbulence) .  On  the  other  hand, 
on  a  large  wing  the  different  intensities  of  the  turbulence  at  different 
points  on  the  span  give  rise  directly  to  a  net  rolling  moment,  which  then 
results  in  rolling  motion.  In  this  section  the  mean-square  value  of  this 
moment  is  calculated. 

At  any  instant  t  the  rolling  moment  L!(t)  can  be  written  in  the 
same  form  as  the  lift  L(t)  in  equation  (4);  however,  according  to  the 
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previously  mentioned  reciprocity  theorem,  the  lift-influence  func¬ 
tion  h(t,y)  is  now  the  lift  distribution  for  an  indicial  roll  with 
unit  helix  angle  at  the  wing  tip.  If  the  assumption  of  invariance  of 
this  distribution  with  time  is  made,  as  for  the  symmetric  case  (see 
e<l*  (5))>  then  the  required  lift-influence  function  can  be  written  as 

h(t,y)  =  i  h'(t)  r'(y) 


where  h!(t)  = 


C^(t)qSb 


2U 


and  where  the  steady-state  lift  distribution 


i  /  \  _  cc 

7  )  =  =^—  now  pertains  to  a  unit  linear  antisymmetric  angle  of  attack. 

Z 

The  correlation  function  for  the  moment  can  then  be  written  as 


t  (  T  )  = 


u  -OO  ^  —00 


'(h)  h'(t2)  VlT+Vts))  dtl  dt2  (22) 


where 


V(Ut)  =  b2  f_^/2  /_b/2  7’(yi)  7'(y2)  y??±f)  ^1 

=  b  f  ^w('lu2'r2+Ti2)  dTl 


(23) 


where,  in  turn,  in  analogy  with  equation  (12), 

(b/2)  -T| 


b/2 


y'(y )  /'(y+n)  dy 


Hence,  the  second  and  third  approaches  indicated  in  the  section  concerned 
with  the  mean-square  lift  (see  eqs.  (15)  and  (l 6))  can  be  used  to  obtain 


26 


NACA  TN  3910 


lr 

and 


where 


and 


dA 


(25) 


so  that  the  mean-square  rolling  moment  can  he  obtained  by  integrating 
this  spectrum. 


A  qualitative  indication  of  the  effect  of  span  on  the  mean-square 
rolling  moment  sensed  by  an  airplane  may  be  obtained  from  the  quantity 

\|r  ?(o),  which  represents  the  integral  of  the  spectrum  CD  ».  For  a 

we  we 

linear  loading,  7  ’  =  6  and 

b/2 


r'(n) 


r’(A) 


6i 


sin 


Ab 

2 


Hence,  using  the  aforementioned  expression  for  the  point  correlation 
function  yields 


q>We-(co) 


-32  +  2P'2  +  (i6(3'2 

(323'  +  6p'3)  K^p')]  +  k'^ J32  -  6p'2  - 
(323'  +  2P'5)  Kx(  p ' )  +  p’3  Ki0(p  *  )j| 


L*w^ 

rtU 


18 


2\ 


+  p'4)  KqCp'  )  + 
I6p'2  KqCp1)  - 


V 


* 
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vhere 


(3*  ==  pj|l  +  kl2 


and 


*we'(°)  =  ^T"  V  “  6  +  e~P(6  +  6P  +  3P2  +  ^  P5) 
0 


This  spectrum  and  its  integral  approach  zero  linearly  as  p  tends  to 
zero. 


Generalized  Aerodynamic  Influence  Functions  in  Unsteady  Flow 

The  aerodynamic  influence  functions  used  in  the  preceding  sections 
define  the  contribution  of  a  given  station  of  a  wing  to  the  total  lift 
and  rolling  moment.  In  the  analysis  of  a  flexible  wing,  and  even  in  the 
calculation  of  certain  properties  of  a  rigid  wing,  generalized  aerodynamic 
influence  functions  are  required,  which  define  the  contribution  of  one 
station  on  the  wing  to  the  lift  at  another  station  and  thus  represent  a 
Green's  function  for  the  unsteady  spanwise  lift  distribution.  No  work 
appears  to  have  been  done  on  such  functions.  For  steady  flow,  apart  from 
some  calculations  for  supersonic  speeds  which  are  based  on  the  subdivision 
of  a  given  wing  into  a  number  of  squares,  the  only  available  results 
appear  to  be  those  given  in  references  9  and  10. 

The  analysis  in  this  section  is  based  on  reference  8  and  consists 
in  a  generalization  of  the  method  presented  therein  to  unsteady  flow. 

This  method  constitutes  an  attempt  to  predict  the  lift  distribution  for 
any  given  twist  on  the  basis  of  knowledge  concerning  a  few  definite 
angle-of-attack  distributions,  and  may  therefore  be  termed  a  function- 
interpolation  method.  For  the  present  purpose,  the  presumably  known  lift 
distributions  are  the  ones  for  uniform  angle  of  attack  in  direct  and 
reverse  flow  in  a  dimensionless  form,  namely,  cc^/cC^;  they  will  be 

referred  to  as  7D(y)  and  yR(y),  respectively.  (The  function  y(y) 
used  previously  is  the  one  now  designated  by  7R(y).)  Also  required  is 
the  coefficient  of  damping  in  roll  C2p .  (The  lift -curve  slope  and  the 

coefficient  of  damping  in  roll  are  the-  same  in  direct  and  reverse  flow 
by  virtue  of  the  reciprocity  theorem,  so  that  no  distinction  will  be 
made . ) 

The  approach  of  reference  3  then  yields  the  following  approximate 
expression  for  the  lift  distribution  due  to  any  angle-of-attack  distribu¬ 
tion: 
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(26) 


(27) 


Values  of  K  may  be  obtained  from  the  information  given  in  references  9 
and  10.  As  the  aspect  ratio  tends  to  zero,  K  approaches  1/2,  whereas 
for  aspect  ratios  approaching  infinity,  K  tends  to  1.  The  following 
relations  can  be  obtained  from  elementary  definitions  and  from  the 
aforementioned  reciprocity  theorem: 


With  the  aid  of  these  relations  the  lift  distributions  given  by  equa¬ 
tion  (26)  may  read.ily  be  seen  to  have  the  correct  lift  and  rolling  moment 
for  all  angle -of-attack  distributions,  and  to  reduce  to  the  exact  lift 
distributions  for  angle -of-attack  distributions  which  vary  linearly  along 


MCA  TN  3910 


29 


the  span.  A  lift  distribution  which  possesses  these  properties  could 
readily  be  obtained  by  approximating  any  given  angle -of -attack  distribu¬ 
tion  by  a  suitably  chosen  linear  one.  However,  this  procedure  would 
yield  poorer  approximations  than  the  ones  furnished  by  equation  (26); 
for  a  parabolic  angle-of -attack  distribution  on  a  wing  of  very  small 
aspect  ratio,  for  instance,  this  procedure  would  yield  the  lift 
distribution 


with  a  bending  moment  too  low  by  17  percent  compared  with  the  one  of  the 
exact  lift  distribution, 


whereas  equation  (2 6)  yields 


with  a  bending  moment  8  percent  too  high.  For  very  large  aspect  ratios, 
7d  and  7R  both  tend  to  the  chord  distribution  c/c,  so  that  the  lift 

distributions  given  by  equation  (26)  tend  to  the  correct  limiting  value. 


(28) 


Therefore  the  accuracy  of  the  results  furnished  by  equation  (26)  may  be 
expected  to  increase  as  the  aspect  ratio  increases,  whereas  the  accuracy 
of  the  other  approximation  is  independent  of  aspect  ratio. 


In  this  connection,  it  may  be  mentioned  that  "strip  theory"  consists 
in  using  equation  (28)  for  all  aspect  ratios  and,  hence,  is  not  very 
satisfactory  for  wings  with  medium  and  low  aspect  ratios.  For  instance, 
for  the  case  discussed  in  the  preceding  paragraph,  it  furnishes  a  bending 
moment  which  is  too  low  by  25  percent  for  a  delta  wing,  and  too  high  by 
25  to  100  percent  for  tapered  wings.  (The  two  figures  pertain  to  taper 
ratios  of  1/4  and  1,  respectively. ) 


Substituting  the  value  for  a 
gives  the  following  expression  for 


from  equation  (27)  into  equation  (26) 
ccz/c: 


CC-, 

O 

C 


Or  pW  2  f 

*  vJ_b/2  {My)  - K  ^  +  bK  5 


(y-n)J  rE(n) 


a(ri)  dr) 
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The  expression  in  braces  in  the  integrand  is  the  desired  generalized 
aerodynamic  influence  function  and  will  be  designated  by  G(y,T)).,  so  that 

0(y,T])  =  |yD(y)  -  K  7R(y)  +  bK  6(y-t])j  yR(  tj  )  (29) 

and 


In  the  limiting  cases  of  wings  of  very  low  or  very  high  aspect  ratio, 
and  7^  approach  a  common  value,  say  7.  Thus,  for  wings  of  very 

low  aspect  ratio. 


0(y,il)  =  |[7(y) 


1  b(y-Ti)  7(11) 


(31a) 


where  5  is  the  Dirac  delta  (unit-impulse)  function,  and  for  wings  of 
very  high  aspect  ratio, 


0(y,T})  =  b  S(y-Ti)  j( T])  (31b) 

which  is  the  Green* s  function  associated  with  strip  theory.  Thus,  on 
wings  of  very  high  aspect  ratio  all  the  lift  produced  hy  the  local  angle 
of  attack  at  a  given  station  is  carried  in  the  immediate  vicinity  of 
that  station,  whereas  on  wings  of  very  low  aspect  ratio  much  of  the  lift 
is  carried  elsewhere.  This  tendency  for  a  given  station  to  affect  a 
greater  portion  of  the  wing  as  the  span  decreases  is,  of  course,  to  he 
expected. 

The  preceding  analysis  can  he  applied  to  the  oscillatory  case  at 
a  given  reduced  frequency  k,  as  well  as  to  the  steady  case.  If  the 
assumption  of  invariance  of  normalized  distributions  with  time  or  fre¬ 
quency  is  made,  as  before,  then  7p(y),  7R(y),  and  K  are  independent 

of  frequency,  so  that  equation  (JO)  can  he  written  as 

cr  (k)  qS  rb/2 

*(y;k)  =  - - -  /  G(y,T))  <x(ti)  dT)  (32a) 

b  J  -1/2 

where  G(yyq) ,  defined  as  before,  is  independent  of  k.  By  applying  a 
Fourier  transformation  to  both  sides  of  this  equation  the  following 
relation  is  then  obtained  for  flight  through  continuously  varying  tur¬ 
bulence  (cf.  eq.  (4)),  as  modified  by  the' assumption  stated  in  equa¬ 
tion  (5): 
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p°°  nb/2 

z(y,t)  =  bL^l)  dtl  b  J  j  G(y^)  w(U(t"tl),T1)  dT1 

where  Z(y,t)  is  the  lift  per  unit  span  at  station  y  and  time  t,  and 
where  the  function  h^(t)  is  the  one  used  previously. 

The  correlation  function  for  this  lift  can  then  be  written  as 

,  p00  n00  ,  pb/2  rb/2  fro - o - 


Hy)(T)  ‘  MV  ati  dt2  ^  f'b/lf_h/2  G(y’V  6(y'V  *w(fVvt2)2+(V’ii)2)  aV 

=  -h  r  T  hL  ("M  hL(tg)  dt-L  dtg  V6(u(T+tl“‘t2)'y)  (33a) 

h  «  -oo  ^  -oo 


where  the  function 


^(Ur.y)  =  J"  c^y^p)  G^y,T]2)  dlll  dT]2 

represents  an  effective  correlation  function,  which  when  transformed  into 
the  equivalent  power  spectrum  cp  (co,y)  can  be  used  to  obtain  the  power 

spectrum  for  2(y,t)  and  hence  its  mean-square  value.  Thus 


q,J(y)^  b2\  u 


Cl^qS 


2 

|0(k)|  cpw  (co,y) 


(331) 


Before  the  calculation  of  cpw^(cu,y)  is  discussed,  the  function 
'iVe  ( Ut  j y )  will  he  defined  in  a  somewhat  more  general  form  than  in  the 
preceding  paragraph  in  order  to  anticipate  future  needs,  namely 

tw€(UT^1^2)  =  ^  f  ^/2f  h/2  G(yi^l)  G(y2^2)  tw(ju2T2+(r|2-T)i)2)  dT)L  dr]. 

If  the  assumption  is  now  made  that  /R(y)  and  /^(y)  are  the  same,  as 

is  the  case  for  all  unswept  wings  and  for  wings  of  very  low  or  very  high 
aspect  ratio,  they  can  both  be  identified  with  the  function  y(y)  used 
previously,  so  that 

G(y,n)  =  [(1  -  K)  7(y)  +  bK  6(y-Tj)]7(Tj) 


(3*0 
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V  (UT^y2)  =  i(1  "  k)V(Ut)  +  K(1  "  K)  V*(Ut^i) 


+  K\(|(u2T2+(y2-yl)2)j  VHyl)  7(y2) 

where  \lc,  (Ut)  is  the  correlation  function  calculated  previously  for 
we 

the  averaged  vertical  component  of  turbulence,  and  where 

tWe*(UT,y)  =  rC t])  tw(ju2T2+(y-r1)2)  dr) 

Similarly,  then, 

Cpwe(“'yl'y2)  =  {(1  ‘  K)2cpwe(“)  +  K(1  “  K)]  9we*(“'yl)  + 

V*(“'y2)]  +  ^  ^(“^2^1)}  7(yl)  7(y2)  (35) 

where  cpv  *(to,y)  is  the  Fourier  transform  of  *(UT,y),  so  that 
e  e 

pb/2 

cp  *(m,y)  =  i  /  7(11)  ^(aj,|y-Tj|)  dr)  (36) 

e  D  J  -b/2 

and  thus  represents  an  averaged  form  of  the  spectrum  cpw  introduced 

previously.  (See  eqs.  (l6)  and  ( 19 ) . )  For  uniform  spanwise  loading  and 
the  point  correlation  function  used  previously,  this  function  is 


-  1ST  | 


3k'2[Ki0(a1)  +  Ki0(a2)  - 


U  +  k'2)  L 

al  Ko(al)  "  a2  +  {^l2lCl(al)  +  a2 


where 


al*a2  H  f  +  k 


and  is  shown  in  figure  4  for  several  values  of 

b/2 
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The  more  restricted  form  of  cpv^  required  in  equation  ( 33"b’ )  can 

now  he  obtained  from  the  more  general  form  given  in  equation  (35)  by 
setting  y2  =  y^  =  y  in  the  latter,  to  yield 


-  (1- 


K)2cp- 


(<o)  +  2K(l  -  K)  cpv*Ky)  +  K?q>w(co) 


y2(y)  (37) 


Mean-Square  Bending  and  Pitching  Moments 

When  the  variation  of  the  gust  intensity  along  the  span  is  taken 

into  account,  the  mean-square  lift  1?  is  not  an  adequate  index  of  the 
stresses  in  the  wing,  nor  can  the  stresses  be  obtained  from  the  mean- 

square  lift  distribution  Z^(y).  Instead,  the  mean-square  bending  and 
twisting  moments,  as  well  as  the  .mean-square  vertical  shear,  must  be 
calculated  directly. 

In  a  manner  analogous  to  that  employed  for  the  lift  and  rolling 
moment,  each  of  these  quantities  can  be  expressed  in  terms  of  a  certain 
influence  function  which,  by  means  of  the  reciprocity  theorem,  can  be 
related  to  a  lift  distribution  on  the  wing  in  reverse  flow.  For  instance, 
for  the  root  bending  moment  the  desired  lift  distribution  is  the  one  for 
an  angle  of  attack  which  is  zero  on  one  wing  and  proportional  to  y  on 
the  other,  as  may  be  seen  from  the  fact  that  the  bending  moment  can  be 
expressed  as 

““  'f.h/2  [,<yfl«(y)  f(y>  ^ 

where 

f(y)  =  0  (y  <  0) 

f(y)  3  y  (y  £  0) 

so  that,  according  to  the  reciprocity  theory  for  unswept  wings  in  steady 
or  indicial  flow, 


mb  = 


a(y)  <3y 


Similarly,  for  the  root  shear  the  required  lift  distribution  corresponds 
to  an  angle  of  attack  which  is  zero  on  one  wing  and  uniform  on  the  other. 
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Such  lift  distributions  can  be  calculated  readily.  For  instance, 
for  subsonic  flow  and  unswept  wings  they  can  be  obtained  from  those  given 
in  reference  10 j  the  lift  distribution  for  the  root  bending  moment  is 
one-half  the  sum  of  a  linear  symmetric  and  linear  antisymmetric  lift 
distribution,  and  the  lift  distribution  for  the  root  shear  is  one-half 
the  sum  of  a  lift  distribution  due  to  a  uniform  angle  of  attack  and  a 
lift  distribution  due  to  deflection  of  a  full-span  aileron. 


However,  in  some  cases  such  calculations  may  be  time-consuming, 
and  an  alternative  approach  may  be  desirable.  One  such  approach  consists 
in  synthesizing  the  desired  influence  function  from  the  generalized  lift- 
influence  (Green's)  function  discussed  in  the  preceding  section.  In 
order  to  illustrate  the  use  of  this  approach,  it  is  adopted  in  this 
section. 

The  bending  moment  at  any  station  y  ^0  ^  y  ^  —j  and  at  any 

time  t  can  be  obtained  from  the  lift  distribution  considered  in  the 
preceding  section  as 


rb/2 

MB(y,t)  =  /  (y1  -  y)  l(y'  ,t)  dy' 

Jy 


p00  ^  b  /2  \ 

=  ^J  (y'  -  y)  <3 y'J  hL(tl)dtigJ  ,  G(y',Tl)  w(u(t-tL)r)) 


dr] 


b/2 


■lP  p  u  /  c.  ^ 

=  bj  MM  dtiJ  j  S(y,n)  M(t_t iM)  dTi 


where  the  influence  function  for  the  bending  moment  is 


(38) 


b/2 


M(y,T))  =  i  J  (y ’  -  y)  G(y ’ ^ "n)  <iy' 


(39a) 


so  that,  upon  introducing  the  previously  used  function  for  G(y,r]) 

M(y,n)  =  £(1  -  K)  ML(y)  +  K  Mg(y,T])J  7(11)  (39b) 
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where 


M2(y,T))  =  n  -  y  (n  >  y) 
M2(y,ii)  =0  (n  =  y) 


and,  hence ,  for  the  root  tending  moment, 

i  rb/2 

M]_(0)  =  y'  /(y* )  dy' 


M2(0,t])  =  tj  (n  >  0) 

M2(0,ri)  =0  (n  <  0) 


Hence,  the  correlation  function  for  the  root  tending  moment  is 

p  CO  p  00  p h/2  ph/2  f 

*»B(T)  ‘H  MM  MM  dtl  dt2  ^J.h/J.i/2  \{1  -  K)  V<^> 

K(1  -  K)  M1(0)[M2(0,Tll)  +  M2(0;r)2j]  + 

K2M2(0^l)  ^(O^jy^i)  7( n2)  ^w(ju2 (T+tl-t2 )2+  (Tl2-r>l)2)  dTll  dTl 

(40) 


The  approach  used  in  the  preceding  sections  can  now  te  used  to  obtain  the 
power  spectrum  of  the  root  tending  moment  by  evaluating  the  inner  pair'  of 
integrals  of  equation  (^0),  taking  the  Fourier  transform  of  the  result 
with  respect  to  t,  and  multiplying  the  power  spectrum  obtained  in  this 


manner  ty 


U  , 


1 000 1 


indicated  for  the  lift. 


or  ty  using  one  of  the  alternative  approaches 


36 


NACA  TN  3910 


The  spectrum  corresponding  to  the  inner  pair  of  integrals  is,  for 
K  =  1  and  uniform  loading, 


cpw  «(a»)  = 


nUp 


(l  +  k'2) 


(-64+  2(3  1 2)  +8p,2K0(^-)  +  (32P’+  P'5) 


k,2jp'3Ki0(^-)  +  (64- 6p'2)  -8p,2K0(^  -  (32P'  -  p'5)  K±(^-)] 


2" 

where  3'  =  +  >  as  ^e^ore. 


The  integral  of  this  spectrum  is 


= 


w1 


3 

3P 


(-24  +  +  e-^2  (24  +  12 p  +  3P2) 


Examination  of  \|l.  »i(o)  indicates  that,  although  the  mean-square  lending 

moment  tends  to  decrease  as  the  span  increases,  it  decreases  less  rapidly 

than  the  mean-square  lift,  with  the  result  that  the  effective  lateral 

center  of  pressure  moves  outboard.  Quantitatively,  the  distance  from  the 

plane  of  symmetry  to  the  effective  lateral  center  of  pressure  can  be 

defined  as  the  square  root  of  the  ratio  of  the  mean-square  root  bending 

moment  to  the  mean  square  of  the  lift  on  one  wing,  that  is,  of  the  root 

shear.  Although  these  mean  squares  have  not  been  calculated,  the  square 

root  of  the  ratio  of  u(0)  to  the  corresponding  value  for  the  root 

we 

shear  increases  by  15.5  percent  as  3  increases  from  0  to  infinity,  with 
much  of  the  increase  realized  at  fairly  small  values  of  3 • 


For  a  swept  wing  the  variation  of  the  gust  intensity  along  the  span 
results  in  a  pitching  moment  which  must  be  taken  into  account  in  calcu¬ 
lations  of  the  dynamic  response  of  the  airplane  to  continuous  turbulence. 
This  pitching  moment  can  be  obtained  in  substantially  the  same  manner 
as  the  bending  moment.  Thus,  if  y  is  the  station  of  the  mean  aero¬ 
dynamic  chord, 

rb/2 

M(t)  =  tan  A  /  (y  -  |y|)  Z(y,t)  dy 
J  -b/2 


b/2 


H  P  p  U/ 

=  tanA  ij  MM  dtlJ  M5(tj)  w(u(t-tL), t))  drj 


(41) 
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where 


ph/ 2 

1*5(11)  =  J  /  (y  -  |y|)  G(y,Ti)  <*y 

^  13  J  -b/2 

so  that,  with  the  previously  used  approximation  to  the  Green's  function, 

M3(tj)  =  |(1  -  K)[y  -  2Ml(0)]  +  K(y  -  h  |)  j  7(tj) 

The  correlation  function,  spectrum,  and  mean-square  value  of  the  pitching 
.moment  can  then  he  obtained  in  the  manner  used  in  the  preceding  sections. 


Wing-Tail  Correlation  Effects 

The  tail  strikes  a  given  gust  some  time  after  the  wing  does;  as  a 
result,  a  pitching  moment  arises  which  does  not  exist  in  steady  or  quasi¬ 
steady  motion,  nor  if  the  airplane  is  very  small,  because  then  the  time 
lag  is  insignificant.  This  pitching  moment  can,  for  the  purpose  of  the 
present  paper,  be  analyzed  either  in  terms  of  the  correlation  between  the 
gusts  at  the  wing  and  those  at  the  tail  or,  if  a  time-lag  term  is  included 
in  the  indicial-response  function  of  the  tail,  in  terms  of  the  correlation 
between  the  wing  and  tail  response  functions.  The  first  point  of  view 
serves  to  exhibit  the  effect  under  consideration  more  clearly  and  is 
adopted  first;  the  second  is  more  convenient  and  is  adopted  in  the  subse¬ 
quent  parts  of  this  report. 

In  the  somewhat  artificial  case  of  a  small  wing  and  tail  separated  by 
a  relatively  large  distance,  only  the  distribution  of  turbulence  along  a 
line  (the  flight  path)  rather  than  in  a  portion  of  a  plane  is  needed. 

The  pitching  moment  due  to  the  vertical  component  of  atmospheric  turbulence 
can  then  be  written  as 


where  x-j.  is  the  tail  length,  and  where  h^t)  and  h^-(t)  are  the 

pitching -moment  responses  to  indicial  gusts  hitting  the  wing  and  tail, 
respectively,  at  t  =  0;  both  may  include  unsteady-lift  effects,  and, 
if  downwash  effects  are  to  be  considered,  h^t)  should  include  the 

contribution  to  the  pitching  moment  of  the  tail  lift  caused  by  the  down- 
wash  at  the  tail  associated  with  the  lift  on  the  wing  which  results  from 
the  indicial  gust. 
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The  spectrum  of  this  moment  can  then  he  written  as 


|h*(o))  +  Ht(o)) 


<PW(«>) 


where  the  symbols  R  j  j  and  *  designate,  respectively,  the  real  part 

and  the  complex  conjugate  of  a  complex  number.  Obviously,  when  x^ 

approaches  0  the  second  term  in  the  bracket  in  this  equation  vanishes, 
so  that  the  first  term  represents  the  perfect-correlation  effect,  and 
the  second  represents  the  correction  for  imperfect  correlation. 


In  order  to  furnish  an  estimate  of  the  magnitude  of  the  effects 
under  consideration,  some  calculations  have  been  made  on  the  basis  of 
the  assumption  that  a  real  coefficient  \i  exists  such  that 

Hif(co)  =  P  H^oj) 


which  implies  that  the  attenuation  with  frequency  of  the  contributions 
of  the  wing  lift  and  tail  lift,  respectively,  to  the  pitching  moment  is 
the  same.  The  ratio  p  is  -1  for  neutral  stability,  and  p  >  -1  for 
stable  flight j  it  is  positive  when  the  aerodynamic-center  location 
(tail  off)  is  behind  the  center  of  gravity,  so  that  positive  values  of 
p  are  not  likely  to  be  incurred  with  normal  configurations  and  flight 
conditions . 


For  this  case, 


i  ,  .  1 2  n  .2  / 

OiXt\ 

jH^Coo)  1  (l  +  p)  -  2p(l 

-C0B  u  )J 

and 


where 


tail  alone 


u  +  n)2 


|Ht(co)|2Cf>w(o>)  (l  -  cos  fp)  dm 

p  00 

J  |=t»|  cpw(a>)  &D 
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The  function  $  has  been  calculated  for  several  values  of  its  arguments 
by  means  of  the  lift-attenuation  function  given  in  equation  (7)  and  the 
point  spectrum  used  in  the  preceding  sections,  and  is  shown  in  figure  5- 
Also  shown  is  the  ratio  of  the  mean-square  moments  as  a  function  of  the 
factor  p.  which  in  figure  5  Is  referred  to  as  The  effect  of 

imperfect  correlation  is  seen  to  be  very  large  as  the  condition  of  neutral 
stability  is  approached;  the  entire  pitching  moment  is  then  the  result 
os  instantaneous  differences  in  gust  intensities  at  the  wing  and  tail. 

In  general,  however,  the  tail  length  and  the  span  are  of  the  same 
order  of  magnitude,  so  that  an  analysis  of  the  effect  of  imperfect  corre¬ 
lation  between  the  wing  and  tail  must  take  into  account  the  averaging 
effect  of  the  wing  span.  The  pitching  moment  at  any  instant  is  then 


n  00 

J  ht  (tl)  v(xt'rtI(t“tl)'0)  dtl 


(42a) 


Hence,  the  power  spectrum  of  this  moment  is 

2  2 

%(<»)  =  |h»|  <fVe(o>)  +  |Ht(a>)|  cpv(co)  + 

2RjeiUXt  Hfc(a i>)  H**(a>)J  cpWe*(a>,0)  (42b) 

where  cpi_  (a>)  is  the  averaged  spectrum  of  equations  (ll)  and  (13),  and 
we 

cp  *(a),0)  is  the  spectrum  of  equation  (36)  for  y  =  0.  (it  should  be 
noted  that  h^(t)  is  the  response  to  an  indicial  response  which  strikes 
the  tail  at  t  =  0;  if  it  were  the  indicial  response  to  a  gust  which 

jU) 

strikes  the  wing  at  t  =  0,  the  factor  e  ^  in  the  preceding  equation 
would  not  be  required.) 


DYNAMICS  OF  THE  RIGID  AIRPLANE 


In  this  part  of  the  report  the  motions  of  a  rigid  airplane  subjected 
to  atmospheric  disturbances  are  considered.  This  jpart  thus  applies 
directly  to  those  airplanes  which  fly  at  relatively  low  speeds  and  do 
not  experience  any  significant  structural  deformations,  and  also  serves 
as  a  preliminary  to  the  treatment  of  the  flexible  airplane  in  the  later 
parts . 
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The  assumption  is  made  that  the  motions  are  small  enough  to  permit 
the  use  of  linear  approximations  to  the  resulting  aerodynamic  forces  and 
to  permit  the  linear  superposition  of  these  forces.  The  longitudinal 
degrees  of  freedom  (pitching,  vertical,  and  horizontal  motion)  and  the 
lateral  degrees  of  freedom  (yaw,  sideslip,  and  roll)  can  therefore  he 
considered  separately.  The  first  two  sections  of  this  part  are  concerned 
with  the  longitudinal  motions  of  a  small  airplane.  The  material  presented 
here  is  thus  a  generalization  of  the  single-degree -of -freedom  analysis 
presented  in  reference  2.  Mathematically,  the  dynamic  system  is  now 
described  by  a  set  of  simultaneous  ordinary  differential  equations  rather 
than  a  single  one,  but  the  problem  of  calculating  the  pertinent  transfer 
functions  is  still  one  of  simple  algebra.  The  extension  of  these  results 
to  the  large  airplane  is  effected  in  the  third  section,  using  the  tech¬ 
niques  developed  in  the  first  part  of  this  report.  The  lateral  motions 
of  a  large  airplane  are  considered  briefly  in  the  fourth  section,  and  the 
combination  of  the  stresses  due  to  longitudinal  and  lateral  motion  of  a 
large  airplane  is  discussed  in  the  fifth  section. 


Equations  of  Longitudinal  Motion 

The  equations  of  motion  of  an  airplane  can  be  expressed  in  several 
coordinate  systems .  The  system  which  is  generally  the  most  convenient 
one  for  analyses  of  airplane  stability  and  is  generally  referred  to  as 
stability  axes,  consists  of  body-centered  axes  which  are  normal  and 
parallel  to  the  relative  air  velocity  and  rotate  with  it  as  the  airplane 
pitches  or  yaws.  (See  ref.  11,  for  instance.)  The  aerodynamic  forces 
related  to  this  axis  system  can  be  measured  more  readily  in  wind  tunnels 
than  those  related  to  other  axis  systems.  In  view  of  the  very  close 
relation  of  a  stability  analysis  to  the  problem  considered  here,  these 
axes  will  be  used  in  this  part  of  the  paper,  but  in  the  analysis  of  the 
flexible  airplane  in  the  subsequent  parts  of  the  report  space-centered 
axes  will  be  used,  because  they  are  slightly  more  convenient  for  that 
purpose. 

The  airplane  will  be  considered  to  be  in  steady  level  flight  prior 
to  disturbance.  The  motions  studied  will  be  the  deviations  from  their 
mean  values j  for  instance,  the  angle  0  considered  here  will  be  the 
difference  between  the  disturbed  and  the  initial  value  of  the  angle  of 
pitch.  Hence,  the  motions  and  forces  calculated  by  the  method  indicated 
here  must  be  added  to  their  mean  values  to  obtain  the  total  motions  and 
forces . 

Inasmuch  as,  for  the  purpose  of  a  statistical  analysis,  the  dynamic 
characteristics  of  the  airplane  are  represented  most  conveniently  by  its 
transfer  functions,  attention  will  be  confined  in  this  section  to  sinus¬ 
oidal  gusts  and  motions. 
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For  this  case  the  linearized  equations  of  longitudinal  motion  can 
he  written  as  follows  (see  eqs.  11-193  of  ref.  11,  for  instance): 


ico  -  2^ 

-Zu 

-iaiU 

Zw 

Zu 

> 

-V 

icjo  -  Xu 
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< 
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xw 

Xu 

1 

u 

(43) 
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-  iaMq 
_ 

0 

Mw 

Mu 

The  stability  derivatives  which  appear  in  these  equations  are  defined 
in  table  1  in  terms  of  conventional  aerodynamic  coefficients,  and  the 
numerical  values  are  given  for  the  example  used  in  reference  11.  (The 
value  of  the  mean  chord  is  not  given  in  ref.  11  but  it  is  assumed  herein 
to  be  10  feet  on  the  basis  of  other  information  given  in  ref.  11.) 

In  analyses  of  the  stability  of  a  rigid  airplane  the  quasi-steady 
approximation  to  unsteady-lift  effects  is  usually  made,  in  which,  in 
effect,  the  forces  corresponding  to  a  steady  attitude,  to  constant  dis¬ 
turbance  velocities,  and  to  constant  accelerations  are  considered.  This 
approximation  is  justified  because  the  motions  of  concern  are  generally 
sufficiently  slow.  For  the  same  reason  this  approximation  can  also  be 
made  in  analyzing  the  response  of  an  airplane  to  atmospheric  turbulence. 

However,  in  this  problem  another  type  of  unsteady-lift  effect  occurs, 
namely,  that  related  to  the  forces  directly  attributable  to  the  turbu¬ 
lence.  This  effect  is  here  taken  into  account  by  multiplying  the  quasi - 
steady  values  of  the  forces  due  to  gusts  on  the  right  side  of  equa¬ 
tion  (ll3)  by  the  attenuation  function  0(k).  This  procedure  implies  the 
assumption  that  the  airplane  is  small  relative  to  the  scale  of  turbulence, 
inasmuch  as  no  averaging  effects  have  been  taken  into  account;  these 
effects  will  be  discussed  presently.  Also,  this  attenuation  function  is 
strictly  applicable  only  to  the  normal  forces. 


The  unsteady  effects  on  the  drag  are  not  known  because  of  the  rela¬ 
tively  complicated  nature  of  the  mechanism  which  gives  rise  to  drag.  If, 
however,  the  assumption  is  made  that  upon  entry  into  a  sharp-edge  gust 
the  drag  rises  linearly  and  attains  its  steady-state  value  in  the  time 

required  to  travel  N  chord  lengths,  the  drag  equivalent  of  |  0(k)| 
is  the  function 


1  -  cos  2Hk 
2N2k2 

which,  for  N  equal  to  about  5  or  6,  agrees  fairly  well  with  |0(k)|2 
in  the  region  of  main  interest  (k  >  l). 
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The  unsteady  moment  is  also  difficult  to  predict  because  of  the 
paucity  of  knowledge  concerning  unsteady  downwash  effects  for  wings  of 
finite  span.  However,  inasmuch  as  the  wing  lift  contributes  part  of 
the  moment  and,  through  the  mechanism  of  downwash,  determines  to  a  large 
extent  the  moment  contributed  by  the  tail,  the  use  of  the  lift  attenu¬ 
ation  function  for  the  moment  appears  reasonable  for  a  first  approxima¬ 
tion,  and  the  use  of  the  same  function  for  the  lift,  drag,  and  moment 
facilitates  the  analysis. 

For  a  more  refined  analysis,  the  time  lag  between  the  instants  at 
which  the  gust  hits  the  wing  and  the  tail  must  be  taken  into  account, 
not  only  for  large  airplanes,  but  even  for  small  airplanes  if  phugoid 
motions  are  important.  (See  ref.  12.)  Also,  the  lag  in  downwash  should 
be  taken  into  account  for  large  airplanes  and  possibly  also  for  small 
airplanes  in  some  cases.  One  way  of  achieving  this  result  is  indicated 
in  the  next  part  of  the  present  paper;  another  is  discussed  in  refer¬ 
ence  12. 

In  equation  (43)  the  unknown  quantities  wp  and  up  are  the  normal 

and  axial  components  of  the  disturbance  velocities  of  the  airplane  rela¬ 
tive  to  the  free  stream.  Inasmuch  as  the  coordinate  axes  rotate  during 
the  motion,  the  time  derivatives  of  these  quantities  do  not  represent 
the  actual  airplane  accelerations,  which  are  required  in  analyses  of  the 
loads  experienced  by  the  airplane  and  the  degree  of  passenger  discomfort. 
If  the  deviations  from  a  mean  flight  path  are  assumed  to  be  small,  the 
accelerations  normal  and  parallel  to  the  chord  or  longitudinal  axis  of 
the  airplane  are  substantially  the  same  as  the  absolute  vertical  and 
horizontal  accelerations  z  and  x,  which  can  be  obtained  from  the 
relations 

z  =  -Wp  +  U(3 
x  =  -Up  -  g0 

The  transfer  functions  for  these  quantities  can  then  be  obtained  by 
introducing  these  relations  into  equation  (43). 

In  studies  of  the  longitudinal  stability  of  airplanes,  equation  (43) 
is  rarely  solved  in  the  form  given  here.  It  is  usually  reduced  to  two 
equations  with  two  unknowns,  either  up  and  0  (the  phugoid  case)  or 

Wp  and  0  (the  short-period  case),  the  short-period  case  being  usually 

the  one  of  primary  interest.  The  part  of  the  turbulent  energy  contained 
in  the  frequency  range  near  the  phugoid  frequency  is  relatively  small, 
so  that  the  phugoid  case  has  no  significance  for  the  analysis  of  loads 
and  accelerations  resulting  from  atmospheric  turbulence.  Hence,  the 
short-period  case,  which  ignores  the  phugoid  oscillations,  furnishes  an 
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excellent  approximation  to  the  loads  and  accelerations  associated  with 
the  longitudinal  motions  of  an  airplane  in  turbulent  air.  However,, 
another  two-degree-of -freedom  case,  the  one  involving  Wp  and  up,  is 
useful  in  certain  studies  of  the  effects  related  to  the  interaction  of 
horizontal  and  vertical  components  of  turbulence. 

Both  of  these  two-degree  cases  can  be  reduced  to  the  single-degree- 
of -freedom  case  involving  only  z  (or  wp) .  For  airplanes  which  have 

a  large  moment  of  inertia  in  pitch  this  simple  case  furnishes  a  good 
approximation.  It  has  been  studied  in  reference  2,  where  substantially 
the  same  approximations  to  the  unsteady-lift  effects  were  made  as  are 
made  here,  except  that  in  reference  2  apparent-mass  effects  were  included 
(These  effects  are  not  included  in  the  stability  derivatives  used  in 
equation  (43)  because  they  are  usually  small  -  less  than  1  percent  of 
the  mass  of  the  airplane  -  and  are  different  for  each  degree  of  freedom. 
However,  if  desired,  the  apparent  mass  pertaining  to  a  given  degree  of 
freedom  can  easily  be  added  to  the  airplane  mass  in  calculating  the 
stability  derivatives.)  However,  this  approximation  is  more  nearly  valid 
for  calculating  peak  loads  persuant  to  an  entry  into  a  sharp-edge  gust 
than  for  calculating  the  response  due  to  random  turbulence.  Consequently, 
in  the  following  sections,  attention  will  be  confined  to  the  short-period 
case,  although  the  analysis  is  equally  applicable  to  the  other  case  and 
easily  extended  to  the  case  of  three  degrees  of  freedom. 


Solution  of  the  Equations  of  Longitudinal  Motion 

Transfer  functions. -  If  the  degree  of  freedom  pertaining  to  x 
(or  up)  is  ignored,  the  solution  of  equation  (43)  can  be  written  as 


J«1 

Hl(a)) 

H§(cjo) 

r  \ 
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►  = 
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$(cd) 

Bg(oj) 

-U 

where  the  transfer  function  h|(cd)  is  defined  by 


H”(<n)  =  0(k) 


.  w  2  ,  _w,  w 

-y  +  Ez^°  +  cz 

— AqDD  +  BqIcd  +  Cq 


(44) 


(45) 


where,  in  turn,  the  coefficients  are  defined  in  terms  of  the  stability 
derivatives  (see  table  l)  by 
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aq  =  1  4  =  -v 

B0  =  -(Zy  +  Mq  +  mt)  b|  =  ^’(UM*  +  Mq) 

C0  =  MqZ^.  -  UMy  4  =  -U(V^  "  MwV) 

The  transfer  function  4(a))  can  be  defined  similarly  in  terms  of  the 
coefficients 

4  =  V  + 

b£  =  -(V^w  -  MwV) 

V 

%  = 0 


In  these  equations  a  distinction  has  been  made  between  the  values  of  Zy 
and  My  which  occur  on  the  right  side  of  equation  (43)  and  are  here 
designated  by  a  prime  mark,  and  those  on  the  left  side  of  that  equation. 
The  primed  derivatives  pertain  to  the  lift  and  moment  directly  due  to 
gusts,  whereas  the  unprimed  derivatives  pertain  to  the  lift  and  moment 
due  to  airplane  motion;  the  reason  for  this  distinction  is  discussed  in 

u  u  u  u  u 

a  later  section.  Furthermore,  the  coefficients  A-^,  B-^,  C4s,  Aq,  Bq, 

and  Cv  are  "the  same  as  the  coefficients  A”,  B”,  •  •  •  except  that 
Zy,  Zy’,  My,  and  My’  are  replaced  by  Z^,  Z^1,  My,  and  Mu'. 

With  these  transfer  functions  the  mean-square  values  of  z  and  & 
can  be  calculated  from  the  spectra  of  w  and  u  by  using  equation  (l), 

provided  that  the  simultaneous  action  of  w  and  u  is  taken  into 

account.  In  order  to  analyze  this  effect  the  vertical  acceleration  z 
will  be  considered,  but  the  analysis  will  be  applicable  to  9  or  any 
other  characteristic  of  the  airplane  which  responds  to  w  and  u. 
Furthermore,  the  transfer  functions  need  not  be  those  considered  in 
equation  (44),  but  can  be  those  calculated  for  the  three -degree-of- 
freedom  system  or  for  a  flexible  airplane. 

Normal-acceleration  spectrum  and  mean-square  value.-  For  the  present 
purpose  the  indicial-response  functions  h-(t)  and  hg(t),  which  are  the 
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Fourier  transforms  of  the  transfer  functions  (see  eq.  (3a) ),  are  more 
convenient.  In  terms  of  these  indicial -response  functions, 

p  00  p  00 

z(t)  =  /  hj(t1)w(u(t-t1))  dtx  +  /  hg(tx)  u(o(t-t1J j  dt-L 

J  —CO  W  -00 

Then,  if  w(t)  and  u(t)  are  stationary  in  a  statistical  sense,  the 
correlation  function  for  2(t)  can  he  written  as 


%(t)  =  f  f  *55  te)  *w(U(T+tI-t2))  +  ^(*1)  ^2)  ^^(T+ti-ta))  + 

O  _00  U  _00  *- 

^ (to)  hg(t2)  ^(T+ti-ta))  +  hI(ti)  hl(t2)  ^(-r+ti-ta))]  atj.  at2  (16) 
where  ^-s  a  cross  correlation  of  w  and  u  defined  by 


^(Ut)  =  w(x+ut)  u(x+U(t+T)) 


Now,  if  the  turbulence  is  isotropic,  the  mutually  perpendicular 
velocity  components  u  and  w  at  points  in  the  XY-plane  are  statis¬ 
tically  independent,  and  their  cross  correlation  is  zero.  Therefore, 
the  two  terms  in  equation  (k6)  involving  ^to(Ut)  vanish,  and  the  power 

spectrum  of  z  is 

=  |i^(cd)  [  ^^.(co)  +  |h£(cd)|2cpu(u>)  (47) 


so  that,  generally  speaking,  the  power  spectrum  of  a  response  which 
depends  on  both  the  horizontal  and  the  vertical  component  of  turbulence 
is  simply  the  sum  of  the  power  spectra  of  the  two  contributions,  provided 
the  turbulence  is  isotropic.  (This  statement  can  be  shown  to  be  true 
even  if  the  distribution  of  the  gusts  over  the  span  is  taken  into  account.) 


For  the  short-period  two-degree -of -freedom  case,  then,  the  contribution 
due  to  w  is 


'-'0  cak  -  (2C0  -  B02)o>2  +  C02 


r 


%(&) 


and  the  contribution  due  to  u  can  be  obtained  from  the  same  expression, 
but  with  the  subscript  and  superscript  w  replaced  by  u.  However,  the 

ratio  of  the  two  contributions  is  in  the  order  of  where  a  is  the 

trim  angle  of  attack  in  radians,  measured  from  the  zero-lift  conditions. 
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Consequently,  except  at  very  high  lift  coefficients ,  such  as  those  used 
in  landing,  the  contribution  due  to  u  is  usually  negligible  compared 
with  the  one  due  to  w  and  is  disregarded  in  subsequent  sections. 

By  means  of  the  lift-attentuation  function  given  in  equation  (7) 
and  the  point  spectrum  used  in  the  first  part  of  the  paper,  the  inte¬ 
gral  for  z^  has  been  evaluated  (by  using  the  technique  of  partial 
fractions  for  the  integrand)  for  the  example  of  reference  11.  (The 
lift-attenuation  function  of  equation  (7)  has  been  used  despite  its 
shortcomings  in  order  to  facilitate  the  analytic  integration  of  the 
spectra. )  The  results  are  shown  in  figure  6,  as  are  the  results  calcu¬ 
lated  similarly  for  the  three-degree-of-freedom  case,  the  other  two- 
degree-of-freedom  case  (horizontal  and  vertical  motion;  referred  to  in 
the  figure  as  the  zero-pitch  approximation),  and  the  single-degree-of- 
freedom  case  (vertical  motion).  An  examination  of  this  figure  indicates 
that,  at  least  for  this  airplane,  inclusion  of  horizontal  motions  does 
not  affect  the  mean-square  normal-acceleration  response  to  any  significant 
extent.  These  calculations  pertain  to  a  "small11  airplane,  and  no  dis¬ 
tinction  has  been  made  between  the  primed  and  unprimed  values  of  Z^ 
and  Mfl-. 

The  preceding  treatment  of  the  short-period  case  has  the  advantage 
of  using  readily  available  information  concerning  the  characteristics  of 
any  given  airplane.  For  the  purpose  of  trend  studies  a  dimensionless 
form  of  the  transfer  functions  is  preferable.  The  preceding  equation 

for  z2  can  be  written  in  dimensionless  form  (the  contribution  of 
horizontal  gusts  again  being  neglected)  as 

£.  /'jjgJf  1 1  rupol2  ^  ^  U  11:' 

g2  W2/  U2*2Jo  k1>_2(ko2_v2)k2+(kc 2+v2)2  L*? 

(48a) 

and,  similarly. 


(2g/c)2 


f)2  *  4  L 


10001 


U  q>wM 


-  2(k02  -  v2)k2  +  (V5  +  v2) 


2  lV 


dk* 


(48b) 


where  k  is  the  mass  parameter 


8m 

Ci^pSc 


K 
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v  is  the  dimensionless  damping  coefficient 

v  =  lQSe  2 
2U  T,  /o 


Tl/2  is  the  time  t0  damp  to  one-'lia]-f’  amplitude,  which  is  given  "by 


+  Mq  + 


so  that 


v  =  i  L  is£  Vi 

A  2%  CI<x  , 


and  kn  is  the  dimensionless  frequency  of  the  short-period  oscillations, 


•with 


so  that 


0>o  =  jMqZvr  -  UM^  -  £(Zw  +  Mq  +  UMfr): 


_2 

1 1  me 


1  me2  ^q.  +  Cmd' 


Thus,  for  this  two-degree-of -freedom  case  the  dimensionless  mean- 
square  responses  are  functions  of  only  two  additional  parameters ,  which 
are  dimensionless  forms  of  the  main  characteristics  of  the  short-period 
case  (the  short-period  frequency  and  the  time  to  damp  to  one-half  ampli¬ 
tude)  ,  "beyond  those  encountered  in  the  single-degree-of -freedom  case, 

_  L* 

namely,  the  mass  parameter  k  and  the  scale  parameter  s  =  _i  - 

Bending -moment  spectrum  and  mean-square  value.-  For  a  small  airplane 
the  instantaneous  "bending  moments  at  various  points  on  the  span  are 
proportional  to  the  instantaneous  normal  acceleration.  For  instance, 
the  root  lending  moment  can  he  written  as 


mb  =  t  y  -  if  y  2 


(49a) 
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where  L^  and  1%  are,  respectively,  the  lift  on  and  the  mass  of  the 

wing,  and  where  y  and  y  are  the  lateral  distances  to  the  center  of 
pressure  of  the  lift  on  one  wing  and  the  center  of  gravity  of  the  mass 
of  one  wing.  However,  is  proportional  to  the  lift  L  on  the  entire 

airplane,  which  in  turn  is  proportional  to  the  normal  acceleration,  so 
that 


and 


(49b) 


Hence,  the  spectrum  and  mean-square  value  of  Mg  are  proportional  to  the 
spectrum  and  mean-square  value  of  z,  respectively,  the  constant  of 
proportionality  being  the  square  of  the  quantity  in  the  brackets  of 
equation  (49b). 


Special  Problems  Related  to  the  Longitudinal  Motion 
of  Large  Airplanes 

Single-degree-of-freedom  case . -  In  the  preceding  sections  the  air¬ 
plane  has  been  assumed  to  be  small  in  the  sense  of  this  report,  and 
neither  the  instantaneous  lateral  variation  of  the  intensity  of  turbulence 
nor  the  difference  between  the  instantaneous  intensities  at  the  wing  and 
the  tail  has  been  taken  into  account.  In  this  section  this  restriction 
is  removed  by  introducing  the  aerodynamic  forces  calculated  in  the  first 
part  of  this  report  into  the  dynamic  analysis  of  the  preceding  sections. 
The  arguments  advanced  in  the  preceding  section  for  ignoring  horizontal 
gusts  and  horizontal  motions  are  equally  valid  for  the  large  airplane; 
therefore  these  gusts  and  motions  will  not  be  considered  here. 

For  the  single-degree-of-freedom  case  involving  only  vertical  motion 
the  required  modification  for  the  normal  acceleration  is  very  simple. 

For  this  case  the  transfer  function  is 

ig(»>  -  -000  7^4-  (50) 

where  ,  attenuated  by  0(k),  is  the  stability  derivative  for  vertical 

gusts  and,  hence,  represents  the  lift  per  unit  gust  intensity.  Therefore, 
if  the  result  for  the  lift  calculated  in  the  first  part  of  this  report 
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is  used,  the  mean-square  normal  acceleration  "becomes 

_ _  2  2 

a2  -  |«k)|2-!^  V(aj)  ^ 

JO  (XT  +  7^d 

(The  function  cp  (o>)  is  defined  in  the  first  part  of  this  report.) 

This  expression  differs  from  the  result  obtained  in  reference  2  only  in 
that  cpw(u>)  is  here  replaced  "by  cpy.  (co) . 


However,  even  for  the  single-degree  case  the  calculation  of  the 
bending  moment  now  becomes  a  considerably  more  complicated  problem, 
because  the  lateral  centers  of  pressure  of  the  lifts  due  to  the  motion 
of  the  airplane  and  directly  due  to  turbulence  no  longer  coincide. 
Equation  (49a)  now  becomes 


Mo  =  -  z  b  4-^7 


y  z  + 


7M( y)  w(u  (t-tjJ  ,yj  dy 


where  is  the  lift-curve  slope  for  the  entire  airplane,  hL(t)  is 

the  response  function  used  in  equation  (j8),  and  7j^(y)  is  the  function 

M(y,T))  used  in  equation  (38),  with  y  =  0.  In  the  first  two  terms  on 
the  right  side  of  equation  (5l),  z  and  z  can  be  related  to  w  by 

means  of  h”  (the  Fourier  transform  of  the  function  E-;  presented 
in  eq.  (50))  and  the  lift  directly  due  to  turbulence  obtained  in  the  first 
part  of  this  paper,  so  that  these  two  terms  can  be  written  in  the  form 

p°°  pb/2 

^M^l)  dtl J  7(7)  v(u(t-t!),y)  dy 

and,  hence ,  equation  (5d)  "becomes 


,(t)  =  if  m  dtl/  [^(M  7u(y)  +  Mh)  7(y)]  w(u(t-t1),y)  dy 
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The  required  influence  function  for  the  bending  moment  is ,  thus, 


^(y^)  =  ^[V1^  7M(y) 


+  hM 


(t)  y(y)] 


(55) 


Although  the  two  terms  of  this  function  are  products  of  a  time-dependent 
and  a  space-dependent  constituent,  as  in  equation  (5);  their  sum  cannot 
he  split  up  in  this  manner.  Hence,  the  techniques  used  in  the  first 
part  of  this  paper  are  not  directly  applicable.  The  general  treatment 
of  problems  involving  influence  functions  for  which  the  assumption  of 
equation  (5)  is  not  valid  will  be  considered  in  the  last  part  of  this 
paper.  However,  in  this  section  a  special  technique  will  be  used  that 
applies  to  cases  for  which  the  influence  function  can  be  expressed  as  a 
sum  of  several  terms  (two  in  this  case),  each  of  which  can  be  expressed 
as  a  product  of  two  functions,  which  depend,  respectively,  on  time  alone 
and  distance  along  the  span  alone.  (See  also  ref.  13- )  This  approach 
is  more  convenient  than  the  general  approach  of  the  last  part  of  this 
report  when  the  number  of  terms  is  two  or,  possibly,  three.  (Although 
in  this  section  only  two  terms  will  be  considered,  the  generalization 
to  three  or  more  terms  is  straightforward. )  When  the  number  of  terms 
is  greater  than  three,  the  general  approach  becomes  more  convenient. 

If  the  value  of  Mg  given  by  equation  (52)  for  t  +  t  is  multi¬ 
plied  by  the  value  for  t  and  the  result  is  averaged,  and  if  the  Fourier 
transform  of  the  resultant  correlation  function  is  then  taken,  the  fol¬ 
lowing  expression  is  obtained  for  the  spectrum  of  the  bending  moment: 


cPmb(^)  =  |hl(co)|  (<»)  +  l\re^  +  28  {hL^  (54) 

where  the  symbols  r|  j  and  *  designate,  respectively,  the  real  part 
and  the  complex  conjugate  of  a  complex  number.  In  this  equation 
and  Hl(co)  are  the  Fourier  transforms  of  hj^(t)  and  hL(t),  cPWe(^) 
is  the  previously  defined  averaged  spectrum  for  the  lift,  and  cpw 

el 

and  (Py.  are  averaged  spectra  obtained  in  a  similar  manner  but  with 
the  value  of  r  given  by  equation  (12)  replaced  by 

(b/2) -r| 

•b/2 


/(y)  7M  (y+O  ay 

s 


(b/2) 


(n)  =  £ J  L7Ms(y)  7Ms(y+1i)  +  ?Ma(y)  ?Ma(y+Tl)J  'V 


(55) 


where  the  subscripts  s  and  a  refer  to  the  symmetric  and  antisymmetric 
parts  of  7^,  respectively. 
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The  contribution  of  the  antisymmetric  part  of  7M  to  ?2  and, 

hence,  to  the  spectrum  for  the  bending  moment  stems  basically  from  the 
asymmetry  of  the  instantaneous  distribution  of  gust  intensity  over  the 
span.  This  asymmetry  gives  rise  to  a  rolling  .moment  (which  was  consid¬ 
ered  in  the  first  part  of  this  report)  and,  hence,  to  rolling  motions, 
which  contribute  additional  bending  moments  due  to  the  aerodynamic  and 
inertia  loads  associated  with  these  motions.  If,  for  the  purpose  of 
calculating  the  bending  moment  due  to  symmetric  flight  through  turbulent 
air  alone,  these  motions  are  disregarded,  then  the  contribution  of  7^ 

to  ?2  should  be  disregarded  as  well.  The  problem  of  combined  symmetric 
and  antisymmetric  motion  will  be  considered  presently. 

As  pointed  out  in  reference  13;  the  mean-square  bending  moment  calcu¬ 
lated  in  this  manner  may  be  smaller  or  greater  than  the  value  calculated 
by  ignoring  spanwise  variations  in  gust  intensity.  If  the  mass  of  the 
airplane  is  almost  entirely  contained  in  the  fuselage,  the  decrease  in 
the  lift  which  results  from  taking  these  variations  into  account  causes 
a  decrease  in  bending  moment.  However,  if  most  of  the  mass  is  in  the 
wing,  the  net  bending  moments  (aerodynamic  less  inertia)  for  a  uniform 
spanwise  gust  are  very  small,  and  the  effect  of  taking  spanwise  variations 
of  gust  intensity  into  account  is  to  increase  the  .mean-square  bending 
moments . 


Two-degree-of -freedom  case.-  For  the  two-degree-of -freedom  (short- 
period)  case  the  analysis  given  for  the  normal  acceleration  in  the  pre¬ 
ceding  section  can  be  extended  as  follows:  As  indicated  in  equation  (44) 
the  transfer  function  for  z  is  now 


■j_  +  (UM^  +  Mqjico  +  UM^ 


m  2  -r>  . 

-u)  +  BqICJQ  +  Cq 


m(-V)  0(k) 


_1_  U(-Zy) 

+  Bq1o>  +  Gq 


lyV  000 


(56) 


y  The  following  analysis  can  be  applied  equally  well  to  0  by  using 
H^(co)  instead  of  H^(to)  .j 


If  the  expressions  inside  the  brackets  of  equation  (56)  are  desig¬ 
nated,  respectively,  by  Hq(o>)  and  Hg^),  and  their  Fourier  transforms 

by  h^(t)  and  h2(t),  then,  as  a  result  of  the  definitions  of 
and  M^1, 


z(t) 


=  /  Mh)  H'-H) 

—CO 


dtf 


1 


h2  tl)  Mt-U)  dtl 


(57) 
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where  L(t)  and  M(t)  are  the  instantaneous  lift  and  pitching  moment 

due  to  the  vertical  component  of  atmospheric  turbulence,  which  have  been 

obtained  in  the  first  part  of  this  report.  The  calculation  of  cp..(to) 

z 

or  of  z  thus  requires  not  only  the  spectra  of  L  and  M,  the  calcu¬ 
lation  of  which  has  been  discussed,  hut  also  the  cross  spectrum  of  L 
and  M,  which  has  to  be  calculated  directly  from  equations  (4)  and  (42a). 
The  result  is 


q>z(o>) 


H1(ul)|2|H(cu)|2cpWe(a))  +  |H2(aj)|2||Hw(co)|\re(a>) 


|Et>)|2cpw(a>)  +  aije^  HfcM  H/(co)  j  ^*(03,0) 
2^  <^(03)  Hg*(a>)  H(a>)  H/U)  }  + 

R  {%(«)  H2*(go)  H(a))  Ht*(a>)J  Cfye*(a3,0) 


(58) 


where  the  first  two  terms  represent  the  contributions  of  the  spectra 
of  L  and  M,  respectively  (see  eqs.  (l4)  and  (42b)),  and  the  third 
represents  the  contribution  of  the  cross  spectrum  of  L  and  M.  For 
the  present  purpose  the  functions  H(cjo),  H^cd),  and  H^(o>)  can  be 
expressed  as 

H(oj)  =  m^)  0(k) 

H^(cd)  =  m  Axa  (-Z^)  0(k) 

Kt(u>)  =  l^yMw  "  m  (-Zw)]  0(k) 

where  Axa  is  the  distance  from  the  aerodynamic  center  (tail  off)  to 

the  airplane  center  of  gravity.  This  definition  of  the  contributions  of 
the  wing  and  tail  to  the  pitching  moment  is  based  on  the  considerations 
that  the  direct  contribution  of  the  wing  can  be  estimated  with  good 
accuracy  and  the  total  pitching  moment  is  likely  to  be  known  from  experi¬ 
ments,  so  that  the  contribution  of  the  tail  (which  includes  the  effect 
of  the  wing  lift  on  the  downwash  at  the  tail)  can  be  determined  as  the 
difference  of  the  two.  The  functions  9We(a>)  and  cpw.e*(a>,0)  have  been 

defined  in  the  first  part  of  this  report.  In  view  of  the  fact  that  the 
function  0(k)  contained  in  some  of  the  terms  of  equation  (58)  always 
appears  in  terms  multiplied  by  others  which  contain  0*(k) ,  only  the 
absolute  square  of  this  function  is  required,  as  before. 


NACA  TN  3910 


53 


Special  Problems  Related  to  the  Lateral  Motion  of 
Large  Airplanes 

The  equations  of  motion  in  the  lateral  degrees  of  freedom  (roll, 
yaw,  sideslip)  have  the  same  form  and  can  he  solved  in  the  same  way  as 
the  equations  for  the  longitudinal  motion.  (See  pp.  111-53  to  III-67 
of  ref.  11.)  Again  it  is  convenient  to  cast  the  problem  in  the  form 
used  in  a  stability  analysis  in  order  to  take  advantage  of  the  results 
of  such  an  analysis.  For  a  small  airplane  it  is  necessary  only  to 
replace  the  terms  due  to  rudder  deflection  by  corresponding  terms  involving 
side  gusts,  namely, 

y45r 

&R 

%R  6r 

where  0'(k)  is  a  suitable  side-force  attenuation  function  for  side 
gusts  and  the  notation  of  reference  11  is  used  for  the  other  terms. 

The  terms  corresponding  to  aileron  deflection  can  be  disregarded  for 
the  small  airplane.  In  the  lateral  degrees  of  freedom  the  small  airplane 
thus  reacts  only  to  side  gusts.  On  the  other  hand,  the  large  airplane 
also  reacts  in  the  lateral  degrees  of  freedom  to  vertical  gusts  through 
the  rolling  moment  calculated  in  the  first  part  of  the  present  paper. 

If  this  rolling  .moment  is  to  be  included,  it  replaces  the  term  5^ 

used  in  reference  11. 

Instead  of  treating  all  three  degrees  of  freedom  simultaneously, 
in  stability  analyses  two  one-degree-of -freedom  cases  are  often  consid¬ 
ered,  namely,  the  one  of  sideslip  alone,  with  angle  of  yaw  equal  and 
opposite  to  angle  of  sideslip  (the  Dutch  roll  case),  and  the  one  of 
rolling  alone. 

The  Dutch  roll  case  may  be  used  for  gust-load  purposes  in  connection 
with  yawing  and  sideslipping  motion  due  to  lateral  gusts,  provided  the 
phase  of  the  motion  is  not  important.  (As  may  be  noted  from  the  preceding 
sections,  the  phase  of  a  transfer  function  is  important  only  in  terms 
involving  cross  spectra. )  Also,  the  Dutch  roll  case  may  prove  useful 
in  calculating  the  vertical-tail  loads  resulting  from  flight  through 
turbulent  air,  particularly  if  the  human  or  automatic  pilot  holds  the 
wings  substantially  level  by  means  of  the  ailerons.  However,  in  general 
it  does  not  appear  to  be  as  satisfactory  an  approximation  as  the  one- 
and  two-degree-of -freedom  approximations  for  longitudinal  motion. 


by  Yvg0’(k) 
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For  rolling  motions  due  to  rolling  moments ,  the  single-degree  case 
of  rolling  alone  appears  to  furnish  a  very  good  approximation.  Although 
the  rolling  .motion  causes  yawing  and  sideslipping  motions,  these  motions 
do  not  appear  to  reflect  on  the  rolling  motion.  Thus,  the  rolling  motion 
which  results  from  the  rolling  moment  can  probably  be  calculated  fairly 
accurately  without  regard  to  the  other  lateral  degrees  of  freedom. 
Furthermore,  within  the  assumption  of  small  motions,  the  stresses  asso¬ 
ciated  with  these  other  lateral  degrees  of  freedom  do  not  generally 
contribute  appreciably  to  those  associated  with  the  longitudinal  degrees 
of  freedom  in  the  parts  of  the  structure  for  which  the  latter  are  critical, 
such  as  the  wing  (although  they  may  be  critical  for  other  parts  of  the 
structure,  such  as  the  vertical  tail).  Therefore  these  degrees  of  freedom 
(yaw  and  sideslip)  will  be  ignored  in  the  treatment  of  the  large  flexible 
airplane  in  the  last  part  of  this  paper.  However,  if  chordwise  bending 
effects  (deformations  parallel  to  the  chord)  are  important,  as  they  may 
be  in  some  cases  at  speeds  close  to  the  flutter  speed,  these  other  lateral 
degrees  of  freedom  may  have  to  be  included  in  the  analysis. 

For  the  large  airplane,  which  responds  in  the  lateral  degrees  of 
freedom  both  to  vertical  gusts  and  to  side  gusts,  the  superposition  of 
the  resulting  responses,  such  as  stresses,  may  be  effected  in  the  way 
indicated  for  the  interaction  of  horizontal  and  vertical  gusts.  If  the 
turbulence  is  isotropic,  the  vertical  and  lateral  gusts  are  statistically 
independent  for  points  in  the  XY-plane,  so  that  the  spectrum  of  a  given 
response  is  equal  to  the  sum  of  the  spectrum  of  that  part  of  the  given 
response  which  is  due  to  vertical  gusts  and  the  spectrum  of  that  part  of 
the  response  which  is  due  to  lateral  gusts. 


Combination  of  the  Results  Obtained  From  Analyses  of  the 
Longitudinal  and  Lateral  Degrees  of  Freedom 

The  instantaneous  wing  stresses  depend  both  on  the  motions  in  the 
longitudinal  degrees  of  freedom  (primarily  vertical  motion  and  pitching) 
and  on  those  in  the  lateral  degrees  of  freedom  (primarily  rolling) .  The 
purpose  of  this  section  is  to  indicate  how  the  stresses  associated  with 
vertical  .motion  and  pitching  can  be  combined  with  those  associated  with 
rolling  due  to  vertical  gusts,  particularly  in  the  case  of  a  large  rigid 
airplane.  (A  small  airplane,  flexible  or  rigid,  does  not  roll  as  a 
result  of  the  action  of  vertical  gusts,  and  for  the  large  flexible  air¬ 
plane  it  is  more  convenient  to  consider  rolling  motion  simultaneously 
with  the  other  motions,  so  that  the  superposition  is  effected  automatically 
in  the  process  of  obtaining  the  required  transfer  functions.)  For  all 
airplanes  the  effect  of  side  gusts  can  then  be  taken  into  account,  if 
isotropy  is  assumed,  by  adding  the  stress  spectra  directly. 


NACA  TN  3910 


55 


In  this  section  the  instantaneous  stress  at  a  given  part  of  the 
ving  will  he  assumed  to  he  proportional  to  the  instantaneous  bending 
.moment  at  that  section  of  the  wing,  so  that  consideration  can  he  confined 
to  this  bending  moment;  also,  the  airplane  will  he  considered  to  he  free 
to  move  in  only  two  degrees  of  freedom,  namely,  vertical  motion  and 
rolling.  The  extension  of  the  following  arguments  to  stresses  which 
depend  on  the  vertical  shear  and  the  torque  as  well,  and  the  inclusion 
of  pitching  as  an  additional  degree  of  freedom,  can  he  effected  readily 
and  will  not  alter  the  conclusion  reached  here. 

When  rolling  motions  of  the  airplane  are  taken  into  account,  the 
bending  moment  due  to  the  aerodynamic  forces  associated  with  rate  of  roll 
and  to  the  inertia  load  associated  with  rolling  acceleration  can  he 
expressed  in  terms  of  the  rolling  moment  by  a  superposition  integral. 

In  turn,  the  rolling  moment  can  he  expressed  in  terms  of  a  superposition 
integral  involving  the  instantaneous  gust  intensities  along  the  span, 
the  influence  function  7*(y)  used  in  the  first  part  of  the  paper,  and 
an  associated  response  function  h^(t).  Hence,  equation  (52)  is  modi¬ 
fied  by  the  inclusion  of  a  third  term  and  becomes 


b/2 


%(*)  =J  dtL  J  [mM  My)  +  MM  ?(y) +  MM  My^^MMM  ^ 

If  the  spectrum  for  Mg  is  now  calculated  in  the  manner  used  pre¬ 
viously,  the  following  result  is  obtained: 
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(59) 


where  H^(od)  is  the  Fourier  transform  of  h^(t),  and  where  cpw  (o>) 

e3 

and  cp^  (go)  can  he  obtained  from  equation  (12)  with  the  following  values 

ek 

of  Tj  and  respectively: 
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(The  function  T^t])  is  four  times  the  function  P 1  ( T] )  considered 

previously  in  connection  with  the  rolling  moment.  This  factor  of  four 
must  he  taken  into  account  in  h^(t).) 


The  spectrum  cp^  may  he  considered  to  consist  of  two  parts:  The 

terms  listed  in  the  first  line  of  equation  (59),  hut  excluding  the  con- 
trihution  of  7^  to  cp^  (see  eq.  (55) )>  represent  the  contribution 

of  the  symmetric  parts  of  the  instantaneous  gust  distributions  or  the 
contribution  associated  with  symmetric  motion;  the  terms  on  the  second 
line  and  the  contribution  of  7^  to  cpv^  represent  the  contribution 

of  the  antisymmetric  parts  of  the  instantaneous  gust  distributions  or 
the  contribution  associated  with  rolling  motion.  Therefore,  the  power 
spectrum  of  the  stress  due  to  gusts  and  combined  symmetric  and  rolling 
.motion  resulting  from  the  gusts  is  the  sum  of  the  two  power  spectra 
(that  for  the  gusts  and  symmetric  motion  alone  and  that  for  the  gusts 
and  rolling  .motion  alone),  provided  the  direct  contribution  of  the  gusts 
is  split  up  into  a  symmetric  and  antisymmetric  part  and  each  is  taken 
into  account  only  once,  in  connection  with  the  appropriate  type  of  motion. 
The  cross-correlation  terms  between  the  symmetric  and  antisymmetric  con¬ 
tributions  to  the  stress  can  be  shown  to  involve  integrals  which  contain 
products  of  symmetric  and  antisymmetric  influence  functions  and,  hence, 
vanish,  so  that  the  cross  correlations  are  zero;  hence,  the  two  parts 
of  the  combined  spectrum  are  statistically  independent  and,  therefore, 
directly  additive. 


DYNAMICS  OF  THE  SMALL  FLEXIBLE  AIRPLANE 


The  purpose  of  this  part  of  the  paper  is  to  consider  the  transfer 
functions  relating  the  stresses  at  various  points  of  a  small  flexible 
airplane  to  the  vertical  gusts  which  cause  them.  The  longitudinal  and 
lateral  degrees  of  freedom  are  still  separable,  and  only  the  longitudinal 
degrees  will  be  considered;  the  lateral  degrees,  which  are  involved  in 
analysis  of  the  response  to  side  gusts  acting  on  the  vertical  tail,  can 
be  analyzed  in  the  same  way.  Therefore,  for  this  case  only  one-half  of 
the  wing  need  be  considered  as  a  result  of  the  symmetry  (or  antisymmetry, 
in  the  case  of  the  lateral  degrees  of  freedom)  of  the  problem. 

The  method  which  is  outlined  in  this  part  consists  in  an  extension 
of  the  numerical-integration  method  of  aeroelastic  analysis  described 
in  reference  Ik  to  sinusoidal  motions  of  the  airplane.  This  extension 
takes  into  account  the  facts  that  the  aerodynamic  forces  now  have  out- 
of -phase  as  well  as  in-phase  parts  and  that  vertical  deflections  must 
now  be  calculated  separately  because  the  structural  deformations  can 
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no  longer  be  characterized  by  angle -of -attack  changes.  Also,  the  "rigid- 
body"  degrees  of  freedom  (vertical  and  pitching  motion  of  the  airplane 
as  a  whole  and  structural  deformations  of  the  tail)  are  now  taken  into 
account. 

The  result  is  a  set  of  linear  algebraic  equations  (which  serve  as 
an  approximation  to  the  ordinary  differential  equations  that  characterize 
the  problem)  for  the  airplane  motions  and  deformations  in  terms  of  the 
applied  aerodynamic  forces  directly  due  to  gusts.  The  desired  transfer 
functions  can  then  be  obtained  from  solutions  of  these  equations  at 
various  frequencies,  and  the  power  spectra  of  the  stresses  are  given  by 
the  product  of  the  absolute  square  of  these  transfer  functions  and  the 
point  power  spectrum  of  the  vertical  component  of  turbulence.  The  same 
equations  may  be  used  to  calculate  with  little  additional  effort  certain 
aeroelastic  effects,  such  as  the  static  aeroelastic  deformations  and  the 
flutter  speed,  which  are  usually  obtained  in  separate  analyses. 


Loads  Applied  to  the  Wing 


The  loads  applied  to  the  wing  stem  from  three  sources:  The  aero¬ 
dynamic  loads  directly  due  to  the  action  of  the  gusts,  the  aerodynamic 
loads  due  to  the  motions  of  the  airplane,  and  the  inertia  loads. 

The  lift  and  pitching  moment  (about  the  elastic  axis)  per  unit  span 
on  a  two-dimensional  airfoil  undergoing  sinusoidal  angle- of -attack  changes 
and  vertical  motions  in  incompressible  flow  are  (see  ref.  15) 


The  terms  multiplied  by  C(k)  are  referred  to  as  the  circulatory  terms 
because  they  are  calculated  from  the  bound  and  shed  vorticity,  and  the 
others  are  referred  to  as  the  potential  terms .  The  potential  terms  are 
in  the  nature  of  additional-apparent-mass  effects,  and  all  those  that 

involve  k2  are  usually  treated  together  with  the  inertia  forces  rather 
than  with  the  aerodynamic  forces.  For  compressible  flow,  however,  the 
forces  are  calculated  in  a  different  manner,  and  the  division  of  the 
forces  into  circulatory  and  potential  parts  then  has  little  meaning. 
Consequently,  in  order  to  facilitate  the  extension  of  this  analysis  to 
compressible  flow,  this  distinction  will  not  be  made  herein. 
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The  aerodynamic  forces  are  therefore  written  as 

n. 


la  =  2rtqc|C-|/k)  a  +  C2(k) 


c/2 


J 


m, 


=  2rfqc2  C5(k)  a  +  Cj^(k)  -5- 


so  that  for  incompressible  flow 

CL(k)  =  (1  +  2e4ik)  C(k)  +  ^  -  e5k2 

C2(k)  =  -ik  C(k)  +  ^ 

C5(k)  =  (l  +  2elfik)e1  C(k)  -  ik  +  +  e^jk2 

C4(k)  =  -ikei  C(k)  -  ^  k2 


(6l) 


In  order  to  calculate  the  lift  at  a  given  point  of  a  wing  of  finite 
span  an  appropriate  Green* s  function  is  required.  An  approximation  to 
this  function  based  on  a  reciprocity  theorem  of  linearized  lifting-surface 
theory,  is  given  in  the  first  part  of  this  paper;  as  used  for  the  compu¬ 
tations  of  that  part,  this  function  implies  the  assumptions  that  the 
spanwise  distribution  of  the  lift  for  oscillations  of  the  wing  as  a  whole 
is  substantially  invariant  with  frequency  and  that  this  distribution  is 
the  same  in  direct  as  in  reverse  flow.  Neither  of  these  assumptions  is 
essential  to  the  analysis  but  both,  and  particularly  the  first,  simplify 
it  considerably.  With  these  assumptions,  the  desired  lift  distribution 
is  then  given  by  expressions  of  the  form  of  equations  (32a)  and  (3^). 

For  the  present  purpose,  however,  a  set  of  aerodynamic  influence 
coefficients  is  required,  rather  than  influence  functions.  Such  a  set 
of  coefficients,  based  on  the  same  ideas,  can  be  obtained  readily  by  the 
techniques  used  in  references  9  and-  10.  The  result  may  be  expressed  as 
follows: 

{Z}a  =  Ciucq|ci(k)[Q]{a}  +  C2(k)  (62a) 

where  the  aerodynamic -influence-coefficient  matrix  M  is  defined  by 

W  =  (1  -  K)p]{l}[ljp]  +  K\y\ 
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where,  in  turn,  {l]>  is  a  unit  column  matrix,  and  LiJ  is  a  row  of 

integrating  coefficients  suitable  for  integrating  a  continuous  function 
for  a  range  of  its  argument  from  0  to  1.  Thus,  for  instance,  if  n 
equidistant  points  on  the  semispan  are  considered,  and  n  is  odd,  then 
according  to  Simpson's  rule 


1  1  4  2  4 

n  -  l\3’  3’  3*  3’ 


Very  little  information  is  available  concerning  the  spanwise  distri¬ 
butions  of  the  pitching  moment  on  wings  of  finite  span  in  unsteady  flow. 
By  means  of  the  reciprocity  theorems  an  appropriate  Green's  function 
could  be  estimated  if  the  lift  distribution  for  wings  with  parabolic 
camber  were  known;  however,  such  lift  distributions  do  not  appear  to  have 
been  calculated  for  wings  of  finite  span  in  unsteady  flow.  In  fact, 
relatively  little  is  known  about  moment  distributions  even  in  steady  flow. 
However,  the  available  information  indicates  that-  the  local  center  of 
pressure  does  not  appear  to  be  very  sensitive  to  the  lift  distribution. 
(See  ref.  9,  for  instance.)  It  will  be  assumed  that  this  is  also  true 
in  unsteady  flow  at  a  given  frequency,  and  that,  furthermore,  these  cen¬ 
ters  of  pressure  are  given  by  two-dimensional  theory.  -With  this  assump¬ 
tion  the  moment  can  be  written  as 

{“t}*  =  +  fy(k)]cj[ci]  j  (62b) 


The  lift  and  moment  distribution  due  to  the  gust  can  be  calculated 
in  a  similar  manner.  In  the  following  analysis,  the  magnitude  of  the 
gust  intensity  is  considered  to  be  unity,  and  the  longitudinal  reference 
point  is  the  intersection  of  the  elastic  axis  and  the  wing  root,  so  that 
the  instantaneous  gust  intensity  at  any  station  y  is 


w  =  e 


-ik!i  tan  A 
c/2 


With  this  function  w, 


{4  - 

R-^MwW 


(63) 


(64) 


where  0(k)  is  the  Sears  function,  as  before. 
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Finally,  the  inertia  loads  are 

=  -nfz  +  me2ca 


The  loads  applied  to  the  wing  can  thus  be  written,  in  summary,  as 


[q]  + 

CLa°  Cx(k)  [Q]  - 

k2  0  tsl 
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k  p  fme2c] 
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(66) 


Loads  Applied  to  the  Tail 


The  loads  applied  to  the  tail  are  similar  in  nature  to  those  applied 
to  the  wing,  but  the  tail  experiences  additional  loads  as  a  result  of 
the  downwash  produced  by  the  lift  on  the  wing.  Again,  little  is  known 
about  the  downwash  in  unsteady  flow,  and  even  in  steady  flow  the  downwash 
cannot  be  predicted  accurately  because  of  boundary— layer  effects  on  the 
fuselage  and  the  wing  root.  Consequently,  even  in  steady-flow  analyses 
experimental  results  are  usually  relied  upon. 

Xn  the  following  analysis,  the  assumption  is  therefore  made  that 
experimental  results  are  available  for  steady  flow,  in  the  form  of  the 
downwash  derivative  de/da.  In  order  to  determine  the  attenuation  of 
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this  value  with  frequency,  the  results  of  the  analysis  of  reference  16 
will  he  used.  These  results  indicate  that  the  time  variation  of  the 
tail  lift  due  to  the  downwash  caused  by  the  wing  lift  which  results  from 
a  unit  jump  in  the  wing  angle  of  attack  can  be  approximated  by  an  imme¬ 
diate  jump  in  the  tail  lift  of  -0.1 6  of  the  steady-state  value  and  another 
jump  to  the  steady-state  value  after  the  time  required  to  travel  the  dis¬ 
tance  from  the  45-percent-chord  point  of  the  wing  to  the  quarter-chord 
point  of  the  tail  plus  another  eighth  of  the  chord  length.  Hence,  for 
sinusoidal  angle-of -attack  changes  the  tail  lift  due  to  downwash  is 


W.,.  ■  -CVA 


8e 

5a 


_ikfel+l£r' 

-0.16  +  1.1 6e  W2  ^ 


(l+2e 41k)  c^-  ik  ^ 

(67) 


where  x-^ 1  is  the  distance  from  the  intersection  of  the  elastic  axis 
and  the  wing  root  (assumed  for  this  purpose  to  be  at  the  ^5  percent  point 
on  the  root  chord)  to  the  aerodynamic  center  of  the  tail.  As  pointed 
out  in  reference  1 6,  this  approximation  is  valid  only  for  k  <  0.35; 
however,  this  range  is  adequate  for  the  present  purpose. 


Similarly,  the  downwash  associated  with  the  wing  lift  due  to  sinus¬ 
oidal  gusts  gives  rise  to  a  tail  lift  which,  within  this  approximation, 
is 
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(68) 


(The  additional  lag  represents  the  time  required  to  travel  the  0.6  root 
semi chords  from  the  45-percent-chord  point  of  the  root,  which  is  the 
reference  point  for  the  gusts,  to  the  75-percent-chord  point  of  the  wing 
root,  which  is  assumed  to  he  the  point  governing  the  lift  at  the  wing 
root,  inasmuch  as  it  is  the  centroid  of  the  influence  function  for  the 
chordwise  pressure  distribution. ) 

The  other  aerodynamic  forces  are  those  due  to  the  motions  of  the 
airplane,  those  due  to  the  tail  deformations,  and  those  directly  due  to 
the  gusts.  On  the  basis  of  the  assumptions  made  in  the  preceding  section, 
these  forces  are 
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where 


,1  +  2ik 
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c5(k)  = 

c6(k)  =  -ik  c(j  k)  +  \  %  k 
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and  cx^.  and  zr  are  the  angle  of  attack  and  vertical  displacement  of 
the  airplane  at  the  wing  root. 


The  inertia  load  on  the  tail  is 


+  Az 


or 


(71) 


Here  the  center  of  gravity  of  the  tail  has  been  assumed  to  coincide  with 
its  aerodynamic  center;  in  order  to  remove  this  assumption  it  is  necessary 
only  to  add  (or  subtract)  the  distance  between  the  two  to  1  in  equa¬ 
tion  (71). 


The  normal  forces  on  the  tail  can  then  be  summarized  as  follows: 

Lt  =  q[Pl(k)  zr  +  F2(k)  +  F3(R)  **  +  Vk)  +  F5OO  wt]  (72) 
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The  pitching  .moments  corresponding  to  these  normal  forces  can  he 
obtained  in  the  maimer  employed  for  the  -wing.  However,  inasmuch  as  the 
tail  chord  is  usually  small  compared  with  the  fuselage  length,  the 
travel  of  the  center  of  pressure  of  the  tail  is  small  compared  with  the 
length  x^!.  Hence,  the  center  of  pressure  will  he  assumed  to  remain 
at  the  aerodynamic  center  of  the  tail,  and  the  pitching  moments  are 
then  -x^’  times  the  corresponding  forces,  so  that 

M.fc  =  -qxp^k)  zr  +  F2(k)  cc*.  +  Fj(k)  Ax  +  F4(k)  Az  +  F?(k)  vt]  (73) 


Wing  and  Tail  Deformations 

The  wing  deformations  may  he  calculated  either  from  structural 
influence  coefficients  or  from  the  bending  and  torsion  stiffnesses  of 
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the  wing  used  in  conjunction  with  simple  beam  theory.  The  latter  approach 
will  be  followed  here,  based  on  the  method  of  reference  12. 


The  bending  and  torsion  moments  on  the  wing  structure  may  be  obtained 
by  integrating  the  applied  loads.  If  numerical  methods  are  employed  to 
perform  these  integrations  the  results  may  be  written  as  follows: 


mbI 


(b/2)s[II] 
cos  A 

-sin  A  -[i] 
2l  j 

_ 

1 

1 - 1 

cos  A  |  [i] 

m. 


T 


(7M 


Similarly,  the  deformations  are 


z  -  z 


a 


r 

- 


r 

(  b  y 

\2  cos  A/ 

&!"] 

M 

i  to] 

|  -  |  tan  A[l"]|^ 

% 

Miji 


bt  7 

W  cos^A 

-®5 

L 

M5  1 

\2/  cos  A 

Mt*ie 

:i] 

-sin  A(|)2[l"] 

mr 


(75) 


where  the  integrating  matrices  (see  ref.  12)  perform  the  following 
operations: 


[lj{f}<>J^  f(i')  d|' 

III]  {f}  -O  fU")  <15"  di' 

[i”]{f}  f(t') 

[lln]|f}  -7>  fU")  <U"  cLg* 
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These  integrating  matrices  may  he  based  on  the  trapezoidal  rule,  Simpson !s 
rule,  or  any  other  numerical  method;  the  intervals  chosen  for  | 

(0  ^  |=l)  need  not  be  of  constant  width  unless  a  specific  rule  demands 
a  uniform  spacing. 

The  structural  deformations  of  the  wing  may  then  be  written  in 
terms  of  the  applied  loads  as 


z  -  zr' 

©  ! 

Jzl 

+. 

t 

0 

1 _ 

M 

a  -  cLp 

*  -  1 

L@i 

K 

>  +  q 

_ig.)  J 

(76) 


where  the  submatrices 


©],  [®]>  [©I  [©],  [©]>  and 


designate,  respectively,  the  four  quadrants  and  two  halves  of  the  matrices, 
obtained  by  postmult  ip  lying  the  square  matrix  of  equation  (75)  by  the 
square  and  the  rectangular  matrix  of  equation  (66),  respectively. 


For  the  purpose  of  the  present  analysis,  which  is  concerned  primarily 
with  the  wing  stresses,  the  tail  deformations  are  treated  by  including 
only  the  vertical  displacement  and  angle -of -attack  change  of  the  tail  as 
a  whole  due  to  the  tail  load.  These  quantities  may  be  obtained  from  a 
static  test  which  consists  in  applying  a  concentrated  normal  load  at  the 
aerodynamic  center  of  the  tail  and  measuring  these  deformations.  They 
may  also  be  obtained  from  a  vibration  test  in  which  the  deflection  curve 
of  the  rear  part  of  the  fuselage  in  the  lowest  vertical-bending  mode  is 
measured;  in  this  case  the  desired  spring  constants  can  be  deduced  from 
the  frequency  relation  of  a  simple  mass  oscillator  in  terns  of  the  meas¬ 
ured  frequency  and  of  the  mass  of  the  empennage  (including  the  part  of 
the  fuselage  which  may  be  considered  to  move  with  the  empennage).  In 
the  absence  of  such  tests  these  constants  may  be  calculated  in  an  anal¬ 
ogous  manner. 


The  tail  deformations  may  then  be  written  as 

Ax  =  -K1Lt 
Az  = 


(77) 


so  that,  also. 


Ax 


Az  = 


(78) 
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Inasmuch  as  these  deformations  are  not  independent  of  each  other,  only 
one  need  he  retained  in  the  analysis.  Therefore,  if  Az  is  eliminated 
hy  means  of  equation  (78) ,  Act  can  then  he  obtained  in  a  form  similar 
to  that  used  for  the  wing  deformations  in  equation  (76): 


Ax  =  -qKL  <  F]_(k) 


zr  +  F2(k)  ct^.  + 


— 

K0 

F3(k) 

-  $  Vk> 

Ax> 

qKx  Fs(k)  wt 


(79) 


Equations  of  Motion 


Equations  (76)  and  (79)  are  equations  of  motion  inasmuch  as  they 
describe  balances  of  aerodynamic,  structural,  and  inertia  forces.  In 
fact,  if  the  airplane  fuselage  were  immobile  (0^  =  zr  =  0),  they  would 

be  sufficient  to  calculate  all  unknown  quantities.  However,  if  the 
fuselage  is  free  to  move,  two  additional  equations  are  required  to  obtain 
the  two  additional  unknown  quantities  oc^  and  zr.  These  additional 
equations  are  those  expressing  the  dynamic  equilibrium  of  the  forces 
on  the  fuselage,  namely, 


-  mf  (zr  -  Ax  =  0 

-  (rp2  +  (AxJ^m-fO^  +  Ax  mfZr  =  0 


(80) 


where  the  wing  lift  and  pitching  moment  can  be  expressed  in  terms  of  the 
lift  and  moment  distributions  Z  and  m  as 


hr  = 

M*  = 


> 


(81) 


or 


=  q 


(82) 
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where 


are  the  rows  obtained  by 


postmultiplying  the  rows  of  equation  (8l)  by  the  square  and  rectangular 
matrices  of  equation  (66).  In  equation  (80)  the  fuselage  lift  and  moment 
have  been  neglected;  they  can  easily  be  expressed  in  terms  of  ctj,  and 
zr  and  included,  if  desired. 


The  equations  for  the  tail  deformation  Ax  and  those  for  the 
overall  normal  force  and  pitching  .moment  can  be  combined  with  equa¬ 
tion  (76)  as  follows:  For  the  sake  of  definiteness  it  will  be  assumed 
that  n  stations  on  the  wing  are  considered,  including  the  one  at  the 
root,  so  that  there  are  2n  +  1  unknown  quantities,  and  that  in  the 
column  matrices  defining  applied  loads,  deformations,  and  so  on,  the 
values  at  the  root  of  the  wing  are  written  at  the  top: 


zroot] 


z 

<  CX  > 

Ax 


Ztip 

}  ^oot  f 


atip 

Ax 


The  first  and  (n  +  l)th  equations  of  the  system  defined  by  equa¬ 
tion  (76)  express  only  the  trivial  fact  that  the  structural  deformation 
at  the  wing  root  is  zero.  They  are  replaced  by  equations  (82)  and  (79), 
which  are  adjoined  to  the  system,  to  yield  the  combined  equation  of 
.motion 


0 

-1  1 

-101 

-1001 


0 


0 

-1  1 

-101 

-100 


0 


0 


(83a) 
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If  the  square  matrix  on  the  right  side  of  this  equation  is  designated 
by  H  ,  the  rectangular  matrix  "by  [b]  ,  and  the  quasi -unit  matrix  on 

the  left  side  of  the  equation  "by  [l*]  ,  the  equation  can  also  he  written 
as 


(83b) 


The  matrices  [©]  to  are  the  same  as  the  matrices  [©] 

of  equation  (76),  except  that  the  first  rows  of  the  latter, 


to 

which  are  all  zero,  are  replaced  as  follows: 


First  row  of  - 

Replacement 

Quantity  added  to  the 
leading  element 

© 

First  half  of  |^(ll)J 

m_p  ^ 

Fx(k)  +  — k2 

id) 

® 

Second  half  of 

Ax  Q 

F2(k)  -  -±— ^  k2 

ti) 

© 

First  half  of  |^(l2) 

-*t  n<*>  -  "© 

id 

© 

Second  half  of  |^(l2)J 

L-  V-f 

-*t  F2(k)  *  “rk2  + 

p(c  ] 

2\2J 

© 

L@J 

0 

© 

L©J 

0 
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Also,  the  elements  of  the  last  rows  and  columns  of  the  matrices  [a] 
and  [b]  are  zero,  except  for  the  following: 


K, 


Al,2»l  *  F3(k)  -  ^  F4<k> 


^n+l,2n+l 

A2n+l,l  = 
A2n+l,n+l  : 

A2n+1, 2n+l 


=  "xt 


Fj(k) 


-Kl  F]_(k) 

=  -Kl  F2(k) 


=  -K- 


F,(k) 


|  Fk<k> 


K 


Bl, n+1  =  F5(k) 
Bn+l,n+l  =  “xt  F5(k) 
B2n+1, n+1  =  -Ki  F5(k) 


Solution  of  the  Equations  of  Motion 


For  the  purpose  of  calculating  the  desired  transfer  functions, 
equation  (83h)  may  be  solved  directly  for  a  given  value  of  q  as  a  set 
of  linear  algebraic  equations  with  coefficients  given  by  the  matrix 

[lr]  -  q[A]  and  with  "knowns"  given  by  the  column  matrix 

^where  [A],  [B]  ,  and  are  functions  of  kj.  The  result  is  a 

column  matrix  of  the  unknown  amplitudes  of  the  motions  of  the  airplane. 
If  this  column  is  calculated  for  several  values  of  k  in  the  range  of 
interest,  these  amplitudes,  considered  as  functions  of  k,  are  transfer 
functions  from  the  gust  to  the  motions. 


This  column  matrix  can  be  substituted  into  equation  (66)  and  the 


resulting  column  matrix 


substituted  into  equation  (74),  to  yield 


a  column  matrix  of  bending  and  twisting  moments  which  again,  considered 
as  a  function  of  k,  represents  transfer  functions  from  the  gust  to  these 
moments.  A  set  of  transfer  functions  for  the  vertical  shear  could  be 
calculated  similarly  from  the  relation 


H  *  [iw  |  H{;} 

The  stress  at  any  point  of  the  structure  can  be  assumed  to  be  given 
by  a  linear  superposition  of  the  bending  moment,  twisting  moment,  and 
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vertical  shear  at  the  given  station,  if  elementary  beam  theory  is  used. 
If  elementary  beam  theory  cannot  be  used  because  of  the  interaction  of 
bending  and  torsion  stresses  or  because  of  shear  lag,  the  stress  at  a 
given  point  can  be  expressed  as  a  linear  superposition  of  moments  and 
shears  at  other  stations  as  well  as  the  given  station.  In  either  case, 
the  transfer  function  for  the  given  stress  is  then  produced  by  the  same 
linear  superposition  of  the  transfer  functions  for  the  corresponding 
moments  and  shears. 

It  may  be  noted  that,  at  zero  frequency,  solution  of  equation  (83b) 
yields  the  static  aeroelastic  deformations  and  thus  permits  the  calcu¬ 
lation  of  the  changes  in  the  lift  distribution  and  the  shift  of  the 
aerodynamic  center  that  result  from  static  aeroelastic  action.  Also, 
inasmuch  as  this  equation  completely  describes  the  dynamic  behavior  of 
the  airplane,  the  speeds  at  which  aeroelastic  instability  phenomena 
occur  can  be  calculated  from  it,  although  such  calculations  are  beyond 
the  scope  of  this  paper.  Suffice  it  to  point  out  that  for  such  a  calcu¬ 
lation  the  degrees  of  freedom  of  the  body  must  be  eliminated  first,  as 
a  result  of  the  way  in  which  the  problem  has  been  set  up.  This  elimin¬ 
ation  can  be  effected  readily  by  considering  the  first  and  (n  4-  l)th 
rows  of  [A],  but  with  AX1,  A-^n+1,  An+1^1,  and  A11+1^n+]_  replaced 

by  0.  If  these  rows  are  premultiplied  by 

A11  Al,n+1  1 

^n+1,1  An+l,n+l 


and  used  as  the  first  and  (n  +  l)th  rows  of  a  matrix  which  is  otherwise 
a  unit  matrix,  and  if  this  resulting  matrix  is  referred  to  as  [ln] , 
then  equation  (83b)  can  be  written  for  this  homogeneous  case  as 


The  products  [l!]  [l"]  and  [a]  [l"]  will  now  have  two  null  rows  and 
columns  each,  which  correspond  to  zr  and  a^.  If  these  rows  and  columns 


are  deleted  and  zr  and  cm,  are  deleted  in  J 


z 

a 

Ax 


>,  yielding  a  column 


^  a  S  ,  the  remaining  matrices  are  nonsingular,  so  that  they  can  be  inverted 
AxJ 

and  the  preceding  equation  can  be  written  as 
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[M  -  <i|p00]] 


z 

a 

Aa 


y  4  0  ? 


where  [e]  is  the  identity  matrix,  and 

[do=)]  =  [tTMj'Twht 


(84) 


The  horizontal  braces  designate  the  fact  that  the  null  rows  and  columns 
have  been  deleted. 

Equation  (8 k)  is  in  the  canonical  form  for  the  calculation  of 
eigenvalues.  If  k  is  set  equal  to  zero  and  the  eigenvalues  of  [d(0 )] 

are  calculated  by  iteration,  expansion  of  the  determinant,  or  any  other 
suitable  method,  the  lowest  real  and  positive  one  represents  the  value 
of  the  dynamic  pressure  at  divergence.  For  swept  wings  the  value  lowest 
in  absolute  magnitude  is  usually  negative  and  is  therefore  of  no  practical 
significance,  although  it  is  often  used  as  an  index  of  the  aeroelastic 
behavior  of  the  airplane. 

This  calculation  can  be  repeated  for  various  positive  values  of  k, 
the  first  few  eigenvalues  being  obtained  for  each.  The  results,  which 
will  generally  be  complex,  can  be  plotted  against  k.  When  any  of  the 
eigenvalues  becomes  purely  real,  it  represents  a  dynamic  pressure  at 
flutter,  and  the  corresponding  value  of  k  represents  the  reduced  fre¬ 
quency  at  flutter.  (This  statement  is  true  only  if  the  structural  damping 
is  zero;  such  damping  effects  can  easily  be  included,  but  the  details  of 
the  process  are  beyond  the  scope  of  this  paper.) 


DYNAMICS  OF  TEE  LARGE  FLEXIBLE  AIRPLANE 


For  the  large  flexible  airplane  the  fundamental  proposition  of 
power  spectral  analysis,  that  the  output  power  spectrum  of  a  system  is 
the  product  of  the  absolute  square  of  the  transfer  function  and  the  input 
power  spectrum,  is  no  longer  valid  if  the  input  is  considered  to  be  the 
gust  intensity  at  a  point.  Nor  can  the  output  power  spectrum  be  expressed 
directly  in  terms  of  an  effective  input  spectrum,  as  in  the  case  of  the 
rigid  airplane,  where  this  simplification  resulted  from  the  assumption 
that  the  indicial-response  function  was  expressible  as  the  product  of  a 
function  of  time  alone  and  a  function  of  distance  along  the  span  alone. 

In  the  first  section  of  this  part  of  the  paper  the  statistical 
problems  involved  in  an  analysis  of  the  response  of  a  large  flexible 
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airplane  are  considered.  The  nature  of  the  generalized  transfer  functions 
required  for  this  purpose  is  described ,  and  the  means  whereby  they  are 
combined  with  the  input  spectrum  are  indicated. 

The  second  section  is  concerned  with  an  extension  of  the  method 
outlined  in  the  preceding  part  to  the  case  of  the  large  airplane. 

Although  fundamentally  the  dynamic  aspects  of  the  problem  are  unchanged, 
and  although  the  longitudinal  and  lateral  degrees  of  freedom  can  still 
be  separated,  a  direct  application  of  the  approaches  outlined  in  the 
preceding  parts  of  the  paper  to  the  large  airplane  requires  consideration 
of  the  entire  wing,  rather  than  only  one  half  of  the  wing.  Little  addi¬ 
tional  computing  time  is  then  required  to  treat  the  lateral  and  longitu¬ 
dinal  degrees  of  freedom  simultaneously,  and  the  necessity  of  combining 
the  results  of  two  separate  analyses  is  obviated.  However,  attention 
can  still  be  confined  to  one  half  of  the  wing  by  using  the  technique 
outlined  in  the  discussion  following  equation  (55)  and  in  the  section 
headed  "Combination  of  the  Results  Obtained  From  Analyses  of  the  Longitu¬ 
dinal  and  Lateral  Degrees  of  Freedom."  Basically  this  technique  consists 
in  splitting  the  influence  functions  of  concern  into  symmetric  and  anti¬ 
symmetric  parts  and  using  one  part  for  an  analysis  involving  the  longitu¬ 
dinal  degrees,  and  the  other  in  an  analysis  involving  the  lateral  degrees. 
If  this  approach  is  adopted,  separate  analysis  of  longitudinal  and  lateral 
degrees  of  freedom  is  still  preferable. 


Extension  of  the  Statistical  Approach 


The  power  spectrum  and,  hence,  the  mean-square  values  of  the 
responses  such  as  the  stress  at  a  given  point  on  the  wing  of  a  large 
flexible  airplane  due  to  flight  through  turbulent  air  can  be  calculated 
in  several  ways.  Perhaps  the  most  direct  of  these  consists  in  using  the 
basic  approach  outlined  in  the  first  part  of  this  paper  and  starting 
with  an  expression  for  the  instantaneous  value  of  the  stress  in  terms 
of  a  suitable  indicial-response  influence  function,  namely, 


ff(t) 


ha(ti>y)  w(u(t-ti)^y)  ^  dti 


(85) 


The  function  hp(t,y)  is,  as  before,  an  indicial-response  influence 

function  or  Green’s  function  for  the  partial  differential  equation.  It 
relates  the  stress  as  a  function  of  space  and  time  coordinates  to  the 
applied  loads,  which  are  also  functions  of  space  and  time  coordinates. 

The  essential  difference  between  the  problem  considered  in  this  section 
and  those  treated  in  the  first  two  parts  of  the  paper  is  that  this  func¬ 
tion  can  no  longer  be  expressed  as  a  product  of  a  function  of  t  alone 
and  y  alone.  In  some  cases  -  in  a  modal  approach,  for  instance  -  it 
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may  "be  expressible  as  a  sum  of  several  such  functions ,  and  then  the 
approach  used  previously  for  the  bending  of  a  large  rigid  airplane  free 
to  move  vertically  may  be  adopted,  as  has  been  done  in  reference  13. 
However,  in  this  section  the  case  is  considered  in  which  even  this 
simplification  cannot  be  made. 

The  correlation  function  for  a  can  be  calculated  directly  from 
equation  (85),  and  for  the  case  of  axisymmetric  turbulence  it  is 


*„(t)  =  /-co/eo/b/s/^  ^2^2)  *w(f2(T+tl-t2)2+(y2-yl)2)  ^1  ^2  dtl  dt2 

Hence,  the  power  spectrum  of  a  can  be  obtained  by  calculating  the 
Fourier  transform  of  this  function  and  is 

pb/2  pb/2  * 

J  b/2J  b/2  ^  M^)  ¥v(“'|y2“yl|)  dyl  ^2  (8T) 

where  cp  (ci},T])  is  the  two-dimensional  spectrum  considered  previously, 

H^(o),y)  is  the  Fourier  transform  of  h^(t,y)  with  respect  to  time, 
and  the  asterisk  designates  a  complex  conjugate,  as  before.  Thus,  the 

function  H^(o),y)  is  a  transfer  function  from  sinusoidal  vertical  gusts 

(of  width  dy)  impinging  on  the  wing  at  a  given  station  y  on  the  wing 
to  the  stress  cr,  or  a  Green  !s  function  for  the  ordinary  differential 
equation  (with  the  quantity  a)  as  a  parameter)  that  relates  the  stress 
amplitude  as  a  function  of  the  space  coordinate  y  to  the  amplitude  of 
the  applied  sinusoidal  gusts. 

The  term  (<*bYp)  Hg(o3,y2)  in  equation  (87)  is  complex;  however, 

the  imaginary  part  can  be  ignored  because  it  contributes  nothing  to  the 
integral  as  a  result  of  the  fact  that  cpv  depends  only  on  |  y2  -  Yp  |  • 

In  a  manner  analogous  to  the  one  employed  in  the  first  part  of  this 
paper,  the  double  integral  in  equation  (87)  can  be  evaluated  by  intro¬ 
ducing  the  function  H^(o),t])  which  takes  the  place  of  the  function  ?(t]) 

used  previously  and  is  defined  by  an- autoconvolution  of  H^(o),r|): 
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where  R<[}'  designates  the  real  part.  Hence, 

p  "b  ^ 

cp  (cu)  =  /  h"(o),ti)  cpw(^il)  dr] 

Oq 


(88) 


Another  approach  consists  in  using  the  double  Fourier  transform 
$w(A)  of  the  axi symmetric  input  correlation  defined  in  the  first  part 

of  this  paper.  In  terms  of  this  function  the  correlation  function  (r) 
can  be  written  as 


’*«'>  ■  HXCCfX  *.(P?)  «: 

=  hi  oo  f  ro  e  1  |^ct(X1^2)|2  H2+*22)  dAl  dX2 


2  dy2  dt^  dt2 


where  the  function 

W  -00  J  -13/2 

represents  the  Fourier  transform  (with  respect  to  y)  of  the  function 
H^(u>,y),  that  is. 


AW 


b/2 


■b/2 


e"lA2y  H^(A1U,y)  dy 


(89) 


Hence, 


With  the  approach  outlined  in  the  next  section,  the  function 

can  le  calculated  either  directly  or  indirectly,  by  first  calcu¬ 
lating  the  transfer  function  from  that  gust  to  the  lift  distribution  and 
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then  the  transfer  function  from  the  lift  distribution  to  the  stress. 
For  the  indirect  method. 


Hj(cd)  G(r],y)  dt| 


where  the  function  H^(co)  G(7},y)  is  the  influence  function  for  the  local 

lift  or  the  Green Ts  function  for  the  aerodynamic  problem  involving  sinus¬ 
oidal  gusts  considered  in  the  first  part  of  this  paper;  the  symbols  r\ 
and  y  in  G(rj,y)  are  interchanged,  however,  so  that  the  function  now 
defines  the  contribution  of  a  gust  at  station  y  to  the  lift  at  sta¬ 
tion  i] .  The  transfer  function  H^(o), rj)  relates  the  (sinusoidal)  stress 

at  the  given  point  to  a  unit  concentrated  (sinusoidal)  normal  force 
acting  at  station  t\. 


With  this  indirect  method,  the  power  spectrum  for  the  stress  at  a 
given  point  can  be  calculated  by  starting  with  the  power  spectrum  for 
the  lift  distribution  calculated  in  the  first  part  of  the  paper.  For 
this  approach  cr(t)  may  be  written  as 


c(t)  = 


Z(t-t1;y)  dy  dtL 


so  that 


p00  p°°  p  b/2  pb/2 

J_cJ_oJ_b/2J-b/2  hff(tl,yi)  +z(T+WyL>y2)  ayL  ay2  ^  at. 


and 


b/2  r b/2  7* 


^  (a»>yi)  Ha(“^2)  ^(“^1^2)  ^1  ^2  (91) 

J  -b/2J  -b/2 


where 


tz(T,yi,y2)  s  l(t,yx)  l{ t+t,y2) 


The  Fourier  transform  of  this  correlation  function  is,  then, 
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(See  the  section  headed  "Generalized  Aerodynamic  Influence  Functions  in 
Unsteady  Flow. ")  Hence, 

I  i2  pV 2  eh/2  *  7 

<Pff(a>)  =  |HZM|  J  J  ^  Hff  (o),yL)  H£(a>,y2)  cpW(_(a>)  dyL  dy2  (95) 


If  cp  (oj)  is  given  hy  equation  (55)j  the  double  integral  can  be 
.  we 

expressed  in  terms  of  single  integrals  as  follows: 
cpa(oo)  =  |rf[(a>)|2|u  -  K)2|g(co)|  2cpWe(ao)  +  2K(l  -  K)R  |g(o>)  G'*(a>)}  + 


f  G"(<b,t))  cpw(ai,n)  dTj 
J  n 


G'V)  =  2  J 


(t/2)-T) 


r  7  l 

R  jH'KyU  R 
,y)|  l|Hj(a),y+ti)| 


K  \  H^^y+n)  }  + 


y(y)  r(y+ii) 


before,  r||  designates  the  real  part,  and  ij  |  designates  the 


imaginary  part. 


Equations  (88),  (90),  and  (93)  thus  represent  three  methods  of 
obtaining  the  power  spectrum  of  the  given  stress.  One  requires  a  trans¬ 
fer  function  from  the  local  gust  intensities  directly  to  the  stress ,  an 
auto convolution  of  this  transfer  function,  and  the  spectrum  of  turbulence 
defined  by  equation  (l6)j  the  second  requires  a  two-dimensional  spectrum 
of  the  turbulence  defined  by  equations  (l7a)  and  (17b)  and  a  Fourier 
transform  of  the  aforementioned  transfer  function  with  respect  to  y; 


NACA  TN  3910 


77 


the  third  utilizes  an  autoconvolution  of  the  transfer  functions  from 
local  concentrated  loads  to  the  stress  and  the  spectrum  for  local  lifts 
due  to  vertical  gusts  calculated  in  the  first  part  of  the  paper.  The 
choice  of  approach  depends  to  some  extent  on  the  information  available, 
but  is  largely  a  matter  of  individual  preference. 


The  analysis  in  this  section  has  been  based  on  the  premise  that 
both  halves  of  the  wing  would  be  treated  simultaneously.  Attention  can 
be  confined  to  one  wing  by  using  only  the  symmetric  parts  of  the  influ¬ 
ence  functions  in  an  analysis  involving  the  longitudinal  degrees  of  freedom, 
and  only  the  antisymmetric  parts  in  an  analysis  involving  the  lateral 
degrees  of  freedom.  The  symmetric  part  of  an  influence  function  for  a 
unit  concentrated  load  or  gust  acting  at  station  y  is  the  response  func¬ 
tion  for  two  loads  or  gusts  of  l/2  unit  intensity  acting  at  stations  y 
and  -y,  respectively.  Similarly,  the  antisymmetric  part  is  the  response 
for  a  load  or  gust  of  l/2  unit  intensity  acting  at  station  y  and  an  equal 
and  opposite  load  or  gust  acting  at  station  -y.  If  this  approach  is 
chosen,  some  of  the  integral  expressions  given  in  this  section  assume 
slightly  simpler  forms.  For  instance,  equation  (89)  can  then  be  written 
for  the  symmetric  part  as 


/\V 


COS  V  <(Hu’y)  ay 


and.  for  the  antisymmetric  part  as 


sin 


Agy  H^(AiU,y)  dy 


Calculation  of  the  Required  Transfer  Functions 

Depending  on  which  of  the  methods  outlined  in  the  preceding  section 
is  used,  one  of  two  types  of  transfer  functions  is  required  -  either  the 
one  from  local  gusts  to  the  stress  of  interest,  or  the  one  from  the  local 
lift  to  that  stress.  Both  of  these  functions  differ  in  several  respects 
from  those  considered  in  the  preceding  part  of  this  paper. 

For  either  type  of  transfer  function  the  structural  deformations  of 
the  wing  under  concentrated  loads  as  well  as  under  distributed  loads 
will  be  required,  so  that  the  numerical-integration  schemes  used  in  the 
preceding  part  have  to  be  modified  to  some  extent.  Also,  it  may  now  be 
preferable  to  consider  the  entire  wing  (rather  than  the  one  semispan 
considered  in  the  case  of  the  small  flexible  airplane) ,  so  that  the 
various  transfer  functions  are  asymmetric.  With  the  degrees  of  freedom 
considered  in  the  analysis  thus  doubled,  it  becomes  preferable  to  include 
one  additional  degree  of  freedom,  roll,  rather  than  to  perform  two  sepa¬ 
rate  analyses  for  symmetric  and  antisymmetric  motions  with,  respectively, 
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one  and  two  less  degrees  of  freedom,  and  then  to  combine  the  results.  How¬ 
ever,  if  the  alternative  approach  of  splitting  up  the  influence  functions 
into  symmetric  and  antisymmetric  parts  for  use  in  two  separate  analyses  is 
adopted,  only  one  semispan  need  be  considered,  and  the  results  presented 
in  this  section  can  then  be  simplified  to  a  large  extent. 

Before  discussing  the  modifications  required  to  extend  the  dynamic 
analysis  outlined  in  the  preceding  part  of  the  paper  to  the  large  air¬ 
plane,  it  might  be  pointed  out  that  chordwise  deformations  (deformations 
parallel  to  the  chord)  will  again  be  ignored.  Again,  they  can  readily 
be  included  by  a  straightforward  extension  of  the  approach  used  here  if 
it  is  felt  that  they  may  be  significant  in  any  given  case.  If  they  are 
included,  however,  yawing  and  possibly  also  sideslipping  motions  can 
probably  no  longer  be  ignored,  because  they  may  give  rise  to  large  forces 
in  the  chordwise  direction.  If  the  entire  wing  is  treated,  these  two 
additional  degrees  of  freedom  can  readily  be  included,  and  all  longitu¬ 
dinal  and  lateral  degrees  of  freedom  are  then  treated  simultaneously; 
if  two  separate  analyses  are  performed  for  the  longitudinal  and  lateral 
degrees  of  freedom,  the  symmetric  and  antisymmetric  parts  of  the  influ¬ 
ence  functions  being  used  and  only  one  semispan  being  treated,  these 
additional  degrees  of  freedom  enter  only  into  the  lateral  analysis. 

The  structural  deformation  due  to  local  (concentrated)  loads  can 
be  obtained  in  several  ways.  If  measured  influence  coefficients  are 
used,  they  pertain  precisely  to  such  loads  and,  in  fact,  must  be  modified 
before  they  can  be  used  for  distributed  loads  (see  ref.  12)  so  that  it 
is  necessary  only  to  use  the  unmodified  coefficients. 

If  the  deformations  are  to  be  calculated  in  a  manner  similar  to  that 
employed  in  the  preceding  part,  the  integrating  matrices  must  be  replaced 
as  follows: 

f[i]  *  |[r] 

where 


■W 

=  1 

(q.  >  p) 

■W 

_  1 

2 

(q  =  p) 

]W 

=  0 

(q.  <  p) 

and 


(I)2 111!  by  It11'] 
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where 


Hpq'  =  yq  -  yP  (a  ^  p) 

Ilpq'  =0  (q  <  p) 


The  factor  of  l/2  for  Ipp'  constitutes  an  approximation  which  implies 
fairing  through  a  discontinuity.  If  this  approximation  is  to  he  avoided, 
the  deflections  due  to  unit  concentrated  loads  (the  structural  influence 
coefficients)  can  he  calculated  directly  from  simple  beam  theory,  in 
which  case  the  limits  of  integration  take  care  of  the  discontinuities . 
Thus,  for  instance,  for  an  unswept  wing,  the  normal  deflection  and  twist 
at  yp  due  to  a  unit  concentrated  load  and  torque,  respectively,  at  y^ 

are 


Zpq 


yq  -  y 

El(y) 


dy  dy' 


zpq 


y 


0  Jo 


yq  -  y 

El(y) 


dy  dy'  +  (yp  -  yq) 


yg  -  y 

El(y) 


dy 


(yP  >  yq) 


“pq 


’it  5JC0  ^  (y*  *  y*>) 

The  concentrated  loads  under  consideration  arise  as  follows:  For 
the  transfer  functions  relating  local  lifts  to  the  desired  stress ,  the 
local  lifts  .may  he  considered  to  he  concentrated  loads  of  unit  magnitude, 
associated  with  concentrated  torques  of  magnitude  e^c.  Equation  (83h) 
can  then  he  written  as  (see  also  eq.  (76)) 


(95) 


where  Id  is  a  diagonal  matrix  of  the  values  of  w  defined  hy  equa¬ 
tion  (63),  and  where  the  matrix  M  represents  either  the  four 

influence-coefficient  matrices  for  z  and  a  due  to  concentrated  loads 
and  torques,  or  the  square  matrix  of  equation  (75)  with  modified  inte¬ 
grating  matrices,  as  discussed  in  the  preceding  paragraphs. 
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It  may  be  noted  that  equation  (95)  now  represents  not  one  set  of 
simultaneous  equations  hut  several,  all  having  the  same  coefficients 
hut  different  sets  of  knowns  (as  defined  hy  the  columns  of  the  matrix 


on  the  right  side)  and,  hence,  different  sets  of  unknowns  j  the  columns 


of  the  matrix 


V 


z 

a 

Ax 


.  This  situation  is  due  to  the  fact  that  the  func¬ 


tions  under  consideration  are,  in  effect,  the  responses  of  the  airplane 
as  a  whole  to  sinusoidally  varying  concentrated  loads  and  are  different 
for  each  location  of  the  applied  load. 


Once  this  equation  has  been  modified  to  take  into  account  the 
overall  body  .motions  and  tail  deflections  (as  explained  in  the  preceding 
part)  as  well  as  the  rolling  motions  (as  explained  in  the  following 


paragraphs),  it  can  be  solved  to  yield  the  unknown  values  of 


From 


these  values  the  bending  and  twisting  moments,  as  well  as  the  vertical 
shears ,  can  be  calculated  and  added  to  those  due  to  the  concentrated 
loads.  When  combined  linearly  as  required  for  the  desired  stress,  these 


moments  and  shears  yield  the  desired  transfer  functions  Hj(u>,y) . 


If  the  transfer  function  directly  from  the  local  gusts  to  the  desired 
stress  is  to  be  determined,  the  response  of  the  airplane  to  the  lift 
distribution  induced  by  a  sinusoidal  gust  of  width  dy  acting  at  sta¬ 
tion  y  must  be  calculated.  This  lift  distribution  is  the  Green's  func¬ 
tion  considered  previously.  If  it  is  represented  by  the  relation 


Hj(a>)  G(y,r)),  with  G(y,rj)  defined  by  the  approximation  given  in  equa¬ 
tion  (34),  the  concentrated  loads  arise  from  the  delta  function  in  that 
expression.  The  right  side  of  equation  (95)  becomes, 


q(l  -  K)  h^(cd)  [r] 


[/(y)  7(7)] 
'[eicJ[/(y )  7(3)] 


R  +  qbK  H^(o>)  [r’J 


that  case. 

' 

|M 

H 

where  [ R  >  represents  the  square  matrix  of  equation  (75),  and  [r  '] 

the  one  discussed  in  connection  with  equation  (95)-  Again,  several  sets 
of  simultaneous  equations  are  implied.  Their  solution  (after  modifica¬ 
tion  for  overall  motions  and  tail  deflections)  yields  values  of  z  and 

a  from  which  the  transfer  functions  H^(o>,y)  can  be  calculated. 
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The  extension  of  the  .method  of  the  preceding  part  to  the  calculation 
of  the  deformations  on  hoth  wings  is  straightforward.  Essentially,  dis¬ 
tributed  lifts  and  torques  now  have  to  be  calculated  for  both  wings  and 
integrated  with  matrices  which  can  be  assembled  from  those  used  for  one 
wing  alone.  No  new  problems  arise  in  this  process,  so  that  it  need  not 
be  discussed  further. 

The  inclusion  of  rolling  motion,  however,  is  not  so  straightforward. 
One  method  consists  in  replacing  all  values  of  z  in  equation  (83a)  or 
its  equivalent  by  z  +  Gy  and  then  reducing  the  columns  involving  this 
quantity  by  the  following  relation  (which  assumes  that  the  new  unknown 
quantity,  the  roll  angle  0,  is  listed  at  the  end  of  the  column): 


* 
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An  additional  equation  must  then  be  joined  to  the  set,  namely,  the 
equation  of  equilibrium  in  roll 


q®  -  MDt®  = 0 


(96) 


where  lx  is  the  inertia  in  roll  of  the  fuselage  and  empennage  alone 
inasmuch  as  the  inertia  effects  of  the  wing  are  included  in  Z(y),  and 
Mp  is  the  coefficient  of  damping  in  roll  for  the  empennage.  For  most 
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cases  both  of  these  contributions  are  negligible.  If  they  are  ignored., 
equation  (9 6)  can  be  written  in  matrix  notation  as 

[nj{i}  .  0 

where  J  is  now  a  matrix  which  serves  to  perform  the  integration 

required  in  equation  (96).  This  condition  can  then  be  adjoined  to  the 
other  equations  of  the  set  in  the  same  manner  as  equations  (82)  were 
adjoined  to  the  set  in  the  preceding  part. 

The  result,  again,  is  a  set  of  simultaneous  equations  for  z,  a, 

Aa,  and  ©  from  the  solution  of  which  the  desired  transfer  functions 
can  be  obtained  as  outlined  in  the  preceding  paragraphs.  Also,  as  before, 
once  the  unknowns  zr,  ar,  and  0  are  eliminated  from  the  set,  the 
divergence  and  flutter  speeds  can  be  calculated  by  conventional  matrix 
operations;  these  speeds  will  then  pertain  to  an  airplane  free  to  move 
vertically  as  well  as  in  pitch  and  roll  and,  hence,  will  include  divergence 
and  flutter  speeds  in  antisymmetric  as  well  as  symmetric  modes. 


DISCUSSION 


Some  Implications  of  the  Assumptions  Concerning  the  Nature 
of  Atmospheric  Turbulence 


The  turbulence  was  assumed  to  be  homogeneous  in  order  to  make  the 
problem  stationary  in  the  statistical  sense  and  thus  permit  the  use  of 
the  mathematical  techniques  developed  for  such  problems.  In  a  practical 
sense,  turbulence  can  be  homogeneous  only  in  a  limited  body  of  air.  The 
assumption  thus  implies  that  the  dimension  of  this  body  of  air  along  the 
flight  path  is  large  compared  with  the  distance  traversed  in  the  reaction 
time  of  the  airplane,  which  in  the  case  of  load  studies  is  of  the  order 
of  the  time  to  damp  to  one-half  amplitude,  but  in  the  case  of  motion 
studies  may  be  much  larger.  Obviously,  the  greater  the  body  of  air,  the 
greater  the  reliability  with  which  the  loads  and  motions  can  be  predicted 
(in  a  statistical  sense)  for  one  run  through  it.  In  general,  turbulence 
at  very  low  altitudes,  which  may  be  influenced  significantly  by  the  con¬ 
figuration  of  the  ground,  and  turbulence  in  thunderstorms  may  not  be 
sufficiently  homogeneous  for  the  purpose  of  this  type  of  analysis,  but 
other  types  of  turbulence  are  likely  to  be  substantially  homogeneous  over 
sufficiently  large  distances. 

Isotropy  was  assumed  in  order  to  permit  the  required  two-dimensional 
correlation  functions  to  be  expressed  simply  in  terms  of  the  one -dimensional 
correlation  functions.  For  sufficiently  short  wave  lengths  all  turbulence 
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is  isotropic,  tut  for  long  wave  lengths  it  can  he  isotropic  only  if  it  is 
homogeneous  (both  in  the  plane  of  the  flight  path  and  P^Pendicuiar  to 
it).  The  condition  of  axisymmetry,  which  is  sufficient  for  most  of  the 
1+s  presented  herein,  is  less  restrictive  than  isotropy  inasmuch  as 
it  does  not  specify  the  variation  of  the  characteristics  of  the  turbulence 
in  the  Vertical  direction.  In  practical  problems,  if  the  turbulence  may 
be  assumed  to  be  homogeneous,  the  conditions  of  axisymmetry  and  i=ut™P7 
are  likely  to  be  satisfied  to  a  sufficient  extent  to  permit  the  nse  ot  ttie 
approach  presented  herein  for  all  but  very  long  wave  lengths .  ® 

length  at  which  it  ceases  to  be  valid  depends  on  the  size  of  the  body  of 
air  under  consideration,  being  larger  for  a  large  body. 

Taylor's  hypothesis  (to  the  effect  that  a  space  displacement  Ax 
alone  the  flight  path  may  be  identified  with  a  time  displacement  r  Ax/U 
in°the^gust  correlation  functions)  implies  that  the  variation  in  gust 
int  ens ity  tha t  prevails  along  the  flight  path  at  any  instant  will  remain 
substantially  the  same  until  the  airplane  has  traversed  the  given  body  of 
air  The  required  correlation  functions  for  atmospheric  turbulence  are 
thus  in  the  nature  of  space  correlation  functions  (rather  than  time  corre- 
latio^f unctions )  and  have  been  considered  as  such.  The  statistical  char¬ 
acteristics  of  the  turbulence  are  then  independent  of  the  speed  at  which 
it  is  traversed.  Clearly,  whether  or  not  this  hypothesis  is  valid  depends 
on  Se  frying  speed  of  the  airplane.  On  the  basis  of  present  knowledge 
no  definite  lower  limiting  speed  can  be  quoted.  Iloweverindications  are 
that  the  hypothesis  is  valid  for  flying  speeds  greater  than  about  100  or 
200  feet  per  second.  The  effect  of  finite  flying  speed  on  the  gust  corre¬ 
lation  function  can  be  expected  to  be  most  pronounced  for  large  distances, 
where  the  correlation  is  weak,  so  that  the  effect  on  the  various  spect 
is  likely  to  be  small  and  to  occur  at  the  longest  wave  lengths,  where,  as 
previously  mentioned,  the  spectrum  is  somewhat  uncertain  for  other  reasras 
as  well;  this  effect  is  thus  more  likely  to  be  significant  for  large  than 
for  small  airplanes . 

The  particular  correlation  function  used  herein  for  the  calculations 
of  the  "averaged"  correlation  functions  and  spectra  has  certain  theoretical 
shortcomings  -  primarily  that  the  associated  spectrum  does  not  decrease 
rapidly  enough  for  very  short  wave  lengths.  However,  it  does  appear  o  e 
adequate  to  Represent  the  available  information  concerning  the  spectra  of 
atmospheric  turbulence  (see  ref.  IT,  for  instance),  because  the  behavior 
at  very  short  wave  lengths  is  relatively  unimportant,  inasmuch  as  airplanes 
cannot  respond  to  them,  and  the  behavior  at  very  long  wave  lengths  is 
usually  in  doubt  by  reason  of  the  nonhomogeneity  of  actual  turbulence.  In 
theTintermediate  range  of  wave  lengths,  this  correlation  function  appears 
to  be  quite  satisfactory . 

The  parameter  L*  (the  integral  scale  of  turbulence)  used  herein  is, 
for  practical  purposes,  a  largely  fictitious  quantity,  inasmuch  as  it  is 
proportional,  to  the  values  of  the  gust  spectrum  for  infinite  vave  lengths. 
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which,  in  view  of  the  uncertainties  in  the  values  of  the  spectra  at  large 
wave  lengths,  have  little  physical  significance.  Therefore,  at  present, 
insufficient  information  is  available  to  give  a  value  for  L*  to  be  used- 
in  connection  with  the  numerical  results  calculated  herein,  although  a 
value  of  1,000  to  2,000  feet  appears  to  be  appropriate.  As  more  informa¬ 
tion  concerning  the  spectrum  of  atmospheric  turbulence  becomes  available, 
more  definite  values  can  be  deduced  by  fitting  an  analytical  expression 
of  the  type  used  here  to  measured  results  in  the  range  of  frequencies  of 
primary  interest,  and  then  using  this  expression  as  a  means  of  obtaining 
a  value  of  L*  by  extrapolation  of  the  measured  results  to  infinite  wave 
lengths  (zero  frequency). 


Considerations  Pertinent  to  the  Application  of  Stationary- 
Random-Process  Techniques  to  Gust -Load  Problems 

The  purpose  of  this  section  is  to  point  out  how,  in  principle,  the 
results  of  analyses  of  the  type  outlined  herein  may  be  used  in  overall 
load  analysis  and,  hence,  in  the  design  of  an  airplane. 

Consideration  is  confined  in  this  paper  primarily  to  the  power  spectra 
of  the  motions  and  stresses  of  interest.  As  pointed  out  in  references  1, 

2,  3,  and  11,  for  instance,  a  great  deal  of  statistical  information  of 
direct  interest  can  be  obtained  from  the  power  spectrum.  For  instance,  if 
the  random  process  of  concern  (say,  the  given  stress  as  a  function  of  time) 
has  a  Gaussian  probability  distribution,  the  expected  number  of  peaks  at 
or  beyond  a  given  level  in  a  given  period  of  time  can  be  calculated  very 
simply  from  the  integral  of  the  spectrum  and  its  second  and  fourth  moments. 

The  results  obtained  in  this  manner  pertain  to  continued  flight  in  a 
given  body  of  turbulent  air.  They  have  to  be  generalized  by  determining 
the  likelihood  of  flying  through  turbulence  of  the  given  characteristics. 
(See  ref.  18.)  The  probability  of  exceeding  a  given  stress  level  during 
the  expected  life  of  the  airplane  while  flying  through  atmospheric  turbu¬ 
lence  can  thus  be  calculated  in  straightforward  fashion.  To  this  proba¬ 
bility  must  then  be  added  the  probability  of  exceeding  this  level  in 
maneuvers,  landings,  and,  possibly,  also  in  turbulence  due  to  thunder¬ 
storms,  because  in  view  of  the  possibly  nonhomogeneous  character  of  turbu¬ 
lence  in  thunderstorms  and  the  possibly  nonlinear  nature  of  the  aerodynamic 
forces  incurred  while  flying  through  them,  the  techniques  used  herein  may 
not  be  applicable  to  flight  through  thunderstorms,  and  a  separate  analysis 
may  have  to  be  performed. 

Although  the  available  information  concerning  atmospheric  turbulence 
is  inadequate  to  permit  of  any  definite  conclusion,  the  results  obtainable 
with  the  approach  outlined  herein  may  turn  out  to  be  most  significant  for 
the  prediction  of  the  low-  and  medium-amplitude  stress  cycles  which  are 
important  to  fatigue  studies;  their  validity  for  or  contribution  to  the 
prediction  of  very  severe  loads  remains  to  be  seen. 
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CONCLUDING  REMARKS 


The  statistical  approach  to  the  problem  of  calculating  the  dynamic 
responses  and  the  stresses  of  an  airplane  subjected  to  continuous  random 
atmospheric  turbulence  has  been  extended  in  several  respects;  basically, 
only  the  assumptions  of  linearity,  that  is,  of  small  motions  and  deforma¬ 
tions,  as  well  as  homogeneity  and  axisymmetry  of  the  turbulence  are 
retained . 

The  first  problem  considered  was  the  effect  of  spanwise  variations 
of  the  instantaneous  turbulent  velocities  on  the  lift  and  moments  due  to 
turbulence.  The  mean-square  lift  has  been  shown  to  be  reduced  consider¬ 
ably  if  the  span  of  the  airplane  is  relatively  large  compared  with  the 
integral  scale  of  turbulence.  The  shape  of  the  spectrum  of  this  lift  is 
affected  relatively  little  by  spanwise  variations  of  gust  intensity, 
except  at  very  hi^ti  frequencies,  if  the  decrease  in  the  effective  mean- 
square  intensity  is  taken  into  account.  The  effect  of  sweep  on  the  mean- 
square  lift  and  its  spectrum  has  been  shown  to  be  small  for  wings  with 
a  given  distance  from  root  to  tip. 

If  the  variation  of  the  instantaneous  velocities  is  taken  into 
account,  the  rolling  moment  to  which  the  airplane  is  subjected  can  be 
calculated.  The  mean-square  rolling  moment  has  been  shown  to  be  propor¬ 
tional  to  the  ratio  of  the  wing  span  to  the  integral  scale  of  turbulence 
for  small  values  of  that  ratio.  Similarly,  expressions  for  the  mean- 
square  values  and  the  power  spectra  of  the  local  lift,  the  bending  moments, 
and  the  pitching  moment  have  been  given.  For  some  of  these  forces  the 
required  aerodynamic  information  cannot  be  calculated  by  existing  methods. 
Therefore,  certain  approximations,  based  on  experience  with  steady  aero¬ 
dynamic  forces  and  available  knowledge  concerning  unsteady  forces,  had  to 
be  made  for  the  aerodynamic  influence  functions  in  unsteady  flow. 

The  next  problem  considered  was  the  dynamic  response  of  a  rigid  air¬ 
plane  to  random  turbulence.  This  problem  had  previously  been  treated 
for  the  case  of  an  airplane  free  to  move  only  in  the  vertical  direction 
and  small  enough  so  that  variation  of  the  turbulent  velocities  along  the 
span  could  be  neglected.  In  the  present  paper  the  response  of  an  airplane 
in  three  longitudinal  degrees  of  freedom  was  considered;  calculations  were 
made  which  suggest  that  the  inclusion  of  deviations  from  the  mean  hori¬ 
zontal  motion  is  superfluous  in  gust-load  calculations.  For  the  remaining 
two  longitudinal  degrees  of  freedom,  the  mean -square  normal  and  angular 
acceleration  have  been  shown  to  be  functions  of  only  two  parameters  other 
than  the  mass  ratio  and  scale  parameter  of  the  single-degree-of -freedom 
case,  namely,  dimensionless  forms  of  the  short-period  frequency  and  of 
the  time  to  damp  to  one-half  amplitude.  An  indication  is  given  of  the 
manner  in  which  the  results  obtained  in  connection  with  the  first  problem 
can  be  used  to  extend  this  dynamic  analysis  to  the  case  in  which  variations 
of  the  turbulent  velocity  along  the  span  have  to  be  taken  into  account. 
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The  last  problem  treated  was  the  dynamic  response  of  a  flexible 
airplane,  including  vertical  motion,  pitch,  and,  when  necessary  (as 
when  spanwise  variations  in  gust  intensity  are  taken  into  account),  roll.  w 

Horizontal  and  lateral  (yawing  and  sideslipping)  motions  were  disregarded 
because  they  do  not  generally  affect  the  wing  stresses  due  to  vertical 
gusts.  A  method  which  represents  an  extension  to  the  dynamic  case  of  a 
numerical -integration  approach  to  the  static  aeroelastic  problem  has  been 
outlined  for  the  analysis  of  the  problem  at  hand.  The  modifications 
required  in  the  basic  statistical  approach  and  in  this  method  of  dynamic 
analysis  in  order  to  treat  the  case  in  which  spanwise  variations  of  the 
gust  intensity  are  important  have  been  discussed. 

Although  most  of  this  analysis  has  been  confined  to  the  vertical 
component  of  turbulence,  it  has  been  shown  that  the  simultaneous  action 
of  longitudinal,  vertical,  and  lateral  gusts  on  the  wing  stresses  (with 
due  allowance  for  the  fact  that  vertical  gusts  affect  both  the  longitudi¬ 
nal  and  the  lateral  motions  of  the  airplane)  can  be  taken  into  account  by 
simply  adding  the  power  spectra  of  the  various  contributions,  provided 
the  turbulence  is  isotropic;  the  cross  correlations  or  spectra  have  been 
shown  to  vanish  either  by  the  symmetry  or  antisymmetry  of  the  influence 
functions  involved  or  as  a  result  of  the  statistical  independence  of 
mutually  perpendicular  velocity  components.  * 

The  approach  presented  herein  thus  furnishes  a  foundation  for  the 
prediction  of  the  statistical  properties  of  the  stress  experience  of  a 
given  airplane  once  the  appropriate  statistical  characteristics  of  the 
atmosphere  have  been  determined. 


Langley  Aeronautical  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Langley  Field,  Va.,  November  5,  1956. 
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TABLE  1.-  DEFINITIONS  AND  NUMERICAL  VALUES  OF 


STABILITY  DERIVATIVES 


Derivative 

Definition 

Numerical  value 
for  the  example 

Zry 

O 

+ 

I 

-1.430 

Zu 

-l?K  +  oi) 

-0.0955 

xw 

i(°I.  -  ° :I>a) 

0.0016 

Xu 

+ c») 

-0.0097 

Mw 

q3c  n 

I  u  % 
iyyu 

-0.0235 

M* 

qSc2  „ 

P 

2IyyU2 

-0.0013 

Mu 

+  ^ 

0 

q.sc2 

2IyyU  m1 

-1.920 

W,  lb .  30,500 

U,  fps .  660 

c,  ft .  10 

Altitude,  ft  .  20,000 


(a)  Unswept  wings. 


(b)  Swept  vinss'  gT  cos  A  =  °-5- 

Figure  1.-  The  mean-square  averaged  vertical  component  of  turbulence 
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Figure  3*-  The  power  spectrum  of  the  averaged  vertical  component  of 

turbulence . 
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Figure  4.-  The  power  spectrum  cpWe*(cD,y) . 
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NACA  -  Langley  Field,  Va. 


Figure  6,-  Effect  of  scale  of  turbulence  on  the  mean-square  acceleration  of  the  example  airplane. 
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